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ABSTRACT 

Context. To obtain an accurate description of broad-band photometric star cluster evolution, certain effects should be 
accounted for. Energy equipartition leads to mass segregation and the preferential loss of low-mass stars, while stellar 
remnants severely influence cluster mass-to-light ratios. Moreover, the stellar initial mass function and cluster metallic- 
ity affect photometry as well. Due to the continuous production of stellar remnants, their impact on cluster photometry 
is strongest for old systems like globular clusters. This, in combination with their low metallicities, evidence for mass 
segregation, and a possibly deviating stellar initial mass function, makes globular clusters interesting test cases for 
cluster models. 

Aims. In this paper we describe cluster models that include the effects of the preferential loss of low-mass stars, stellar 
remnants, choice of initial mass function and metallicity. The photometric evolution of clusters is predicted, and the 
results are applied to Galactic globular clusters. 

Methods. The cluster models presented in this paper represent an analytical description of the evolution of the under- 
lying stellar mass function due to stellar evolution and dynamical cluster dissolution. Stellar remnants are included by 
using initial-remnant mass relations, while cluster photometry is computed from the Padova 1999 isochrones. 
Results. Our study shows that the preferential loss of low-mass stars strongly affects cluster magnitude, colour and 
mass-to-light ratio evolution, as it increases cluster magnitudes and strongly decreases mass-to-light ratios. The effects 
of stellar remnants are prominent in the evolution of cluster mass, magnitude and mass-to-light ratio, while variations of 
the initial mass function induce similar, but smaller changes. Metallicity plays an important role for cluster magnitude, 
colour and mass-to-light ratio evolution. The different effects can be clearly separated with our models. We apply the 
models to the Galactic globular cluster population. Its properties like the magnitude, colour and mass-to-light ratio 
ranges are well reproduced with our models, provided that the preferential loss of low-mass stars and stellar remnants 
are included. We also show that the mass-to-light ratios of clusters of similar ages and metallicities cannot be assumed 
to be constant for all cluster luminosities. Instead, mass-to-light ratio increases with cluster luminosity and mass. 
Conclusions. These models underline the importance of more detailed cluster models when considering cluster photome- 
try. By including the preferential loss of low-mass stars and the presence of stellar remnants, the magnitude, colour and 
mass-to-light ratio ranges of modelled globular clusters are significantly altered. With the analytic framework provided 
in this paper, observed cluster properties can be interpreted in a more complete perspective. 

Key words. Galaxy: globular clusters: general - Galaxy: open clusters and associations: general - galaxies: star clusters 
- galaxies: stellar content 



1. Introduction 

In recent studies, the photometric evolution of star clus- 
ters h as been extensively treated from various approaches 
(e.g., Anders & Fritze-v. Alvcnslcbcn '2003\ Lamcr s et al.l 
[2006; von Hippcl ct al. 2006; Fagiolini ct al. 2007). Because 
cluster photometry is used for a broad range of applica- 
tions, like age-dating galaxies and tracking their forma- 
tion history, it is crucial to obtain an accurate descrip- 
tion of the photometric evolution of clusters. While Simple 
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* The models presented in this paper are publicly avail- 
able in electronic form at the CDS via anonymous 
ftp to http://cdsweb.u-strasbg.fr/ (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ and also 
at http : / / www . astro . uu . nl/^krui j s . 



Ste llar Population (SSP) models (e.g., Leitherer et al.lll999l: 
Br uzual fc Chariot! 120031 : lAnders fc Fritze-v. AlvenslebenI 
120031 : iMarastonI I2005D only consider the changing pho- 
tometric properties due to stellar evolution, other mod- 
els that also use the dynamical input of iV-body simu- 
lations can predict the photometric evoluti on of clusters 
unde r a wider variety of conditions (e.g., Lamers et al.l 
120061 : iFagiohni etall 120071 : iBorch et al.l |2007t). In reality, 
not only stellar evolution but also the dynamical interac- 
tion of a cluster with its environm ent causes it to lose stars 
(e.g.. iBaumgardt fc MakhMl 120031) . This process, which is 
called dissolution, occurs due to internal two-body relax- 
ation and external effects like tidal perturbation, spiral 
arm passages or encouters with Giant Molecular Cloud s 
(e.g.- lBaumgardt fc MakiMll2003t Icieles et al.l[2006Ll2007D . 
It can change the shape of the stellar mass function, 
and will also affect photometric cluster evolution. This 
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is the case as a cluster evolves towards energy equiparti- 
tion, causin g it to preferent ially lose low-mass stars (e.g. , 
Portcgics Z wart et al.l 120011 : [Baumgardt & M akinol l2003t 
Hurlev et al.ll2005D . 



The physical driving force of the preferential loss of low- 
mass bodies is subject to debate. On the one hand, en- 
ergy equipartition between the bodies constituting a clus- 
ter increases the velocities of low-mass objects and thereby 
gives rise to the preferential loss of low-mass stars. On 
the other hand, it has been proposed that mass segrega- 
tion, a phenomenon in which due to energy equipartition 
the more massive stars sink towards the cluster centre and 
low-mass objects move outwards, leads to the same effect 
since bodies in the cluster outskirts are more loosely bound 
than objects in the c luster centre and are thus more eas- 
ily lost (e.g.,lLeon et a l. 2000; Portcgi esZwart et al.ir200H 



iLamers et al 



_ _.] [20061). T his line of reasoning is not compat- 
ible with iKind (l966D . where it is shown that the escape 
rate of stars from a cluster does not vary with radius. In 
that scenario, the preferential loss of low-mass stars and 
mass segregation are both the effects of energy equiparti- 
tion, but do not necessarily sh are any causal connection 
(e.g.. Ide la Fuente Marcos! |20000 . Regardless of its specific 
nature, in the remainder of this paper we consider energy 
equipartition to be the fundamental cause of the prefer- 
ential loss of low-mass stars: either directly, via mass seg- 
regation, or a combination of the two. In any case, mass 
segregation can serve as an indicator for clusters that have 
undergone a strong preferen ti al loss of low-mass stars (e.g. , 
iPortegies Zwart et al.l 120011: iBaumgardt fc Makinol I2003D 
and will therefore be used in that respect. 

There have been many observations of Galactic 
open and globular clusters in w hich evidence of 
mass segregation was found (e.g . , [Anderson &: Kin3 



1996t iHillenbrand k. HartmannI 119 98': 'Zoc cah et all 



Albrow et al.l 



Pasauali et al 



20021: iRicher et al.ll2004 : .Koch et all 
12001 . These clusters can be expected to 



1998; 



20041: 



exhibit non-canonical photometric evolution. Furthermore, 
for a number of clusters overall mass-to-light ratios are 
observed that strongly dev iate fr om the mean value pre- 
sented inlMcLaughlinl (|2000D (e.g.. lBaumgardt et al.ll2003bt 
Ivan de Ven et al. I I2006D . which suggests a range of scenarios 
for photometric cluster evolution. This, in combination with 
the high mass-t o-light ratio in the centre of some globu- 
lar cl usters (e.g.. iPasauali et al]|2004lvan den Bosch et al.l 
|2006[) and the consequent in vocation of intermediate mass 
black holes (IMBHs) (e g., IPortegies Zwart fc McMillanI 
I2OO2I: iGiirkan et all |20 04 l), asks for a cluster model that 
can explain the observed range of mass-to-light ratios. 
While some studies suggest IMBHs to explain the high 
mass- to- light ratio in the centres of globular clusters 
(e.g., iGebhardt eral! 120051 : iNovola et all 120061) . others 
show that these are not required and central concentra- 
tions of stellar remnan ts also provide a solution (e.g., 
IBaumgardt et all l2003allbl : iHurlevI [20071 ). Therefore, it is 
important to investigate to what extent either model can 
be used to explain the observed range of mass-to-light 
ratios, colours and magnitudes. A model describing clus- 
ter mass-to-light ratios may also be able to provide in- 
sight in the connection between globular clusters and ultra- 
compact dwarf galaxies (UCDs), the latter having a differ- 
ent mass-to-light ratio range than globular clusters (e.g., 
Hascgan ct al. 2005; Evstignc eva et al.l2007l : [Reikuba et al.l 
2007; Mieske fc Kroupa,2008i) . 



In this paper we present cluster models that are based 
on stellar isochrones like all SSP models, but analyti- 
cally incorporate dynamical effects on cluster photom- 
etry by following the resul ts from iV-body simulations 
([Baumgardt fc Makinoll20d3l ). The resulting speed and ap- 
plicability to a large parameter space makes it very suitable 
for studying the effect of a range of parameters on cluster 
evolution. We show how photometric properties of clusters 
like their magnitude, colour and overall mass-to-light ratio 
are affected by the preferential loss of low-mass stars, the 
inclusion of stellar remnants, the stellar initial mass func- 
tion (IMF) and metallicity. 

The structure of the paper is as follows. Our cluster 
evolution models are presented in Sect. [2] In Sect. 12.21 it 
is shown how stellar evolution affects the cluster content, 
including the production of stellar remnants. We derive the 
equations describing dynamical effects of cl uster evolution 
on a multi-component powerlaw IMF (e.g.. iKroupal [200l[ ) 
in Sect. 12.31 The effects of the preferential loss of low-mass 
stars and the dynamical loss of remnants are included. In 
that section, we also provide the final set of equations to de- 
scribe cluster evolution with our models. The computation 
of photometric evolution is treated in Sect. |3l The results 
are presented in Sect. HI where also the influences of the 
preferential loss of low-mass stars, stellar remnants, IMF 
and metallicity are investigated, and the results are com- 
pared to previous studies. In SectO the models are applied 
to Galactic globular clusters. Section |S| contains a discus- 
sion of the results, while our conclusions are provided in 
Sect. El 



2. Cluster evolution 

In this section wc describe cluster evolution: first we for- 
mulate the stellar mass function in the cluster, secondly 
we describe cluster mass loss due to stellar evolution, and 
thirdly cluster mass loss due to dissolution. We include the 
effects of the preferential loss of low-mass stars and stellar 
remnants. The models presented here represent the clus- 
ter evolution part of our new cluster population synthesis 
code SPAC]^ which is a fast code to predict observables 
like age, mass and magnitude distributions of clusters for a 
range of galactic conditions. 



2.1. The stellar mass function 

Stars in a cluster are distributed according to a mass func- 
tion, for which we assume a general expression for a multi- 
component powerlaw mass function, which is a function of 
time: 



N,{t,m,)dm, - C(i)77(ms)m3-''('"=Mms, 



(1) 



for m,„in(t) < nis < m^s^x{t), where Ns{t,ms) represents 
the number of stars per M0 at age t, C{t) is the (time- 
dependent) normalisation of the mass function, rjinis) is 
introduced to preserve continuity at the stellar mass where 
the slope (3{ms) changes, and Wg denotes stellar mass in 
M0. Setting t — gives the initial mass function (IMF): 
Ns{0,ms). For a Salpeter IMF, we have constant values 



^ This is an acronym for Stellar Population Age Computing 
Environment. 
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ri{ms) — 1 and /3(ms) — 2.35 with mmin.i ~ 0.1 Mq, while 
a Kroupa IMF is described by 



/3(ms) 



Pi = 1.3 for 0.08 Mq < < 0.50 Mq, 
/32 = 2.3 for 0.50 Mq < m^, 



(2) 



with initial minumum stellar mass m^ia^i — 0.08 Mq and 
r/(m, > 0.50 Mq) = 2r]{ms < 0.50 Mq). 

After cluster formation, clusters lose mass due to stellar 
evolution and dissolutior0. This section provides a descrip- 
tion of how both mechanisms affect the stellar mass func- 
tion and cluster content. Clusters with and without the 
preferential loss of low-mass stars are treated. 

2.2. Stellar evolution 

The total mass loss due to stellar evolution and dissolution 
can be written as 



dt 



dM^f 
dt 



dM^f 
dt 



(3) 



dis 



In this section we describe the mass loss due to stellar evo- 
lution. We separate cluster mass into its two constituents: 
the mass in stars, or the luminous mass M^J™(i), and the 
mass in stellar remnants M^^(t). The total mass Af^°*(t) is 
then obtained from 



youngest isochrone {\ogt = 6.6). After assuming an IMF, 
i.e., iVs(0, rus), we can write the remaining mass fraction of 
luminous mass as 



.it) 



msiVs(0,r7is)dms 



cl,i 



(8) 



/xj,™(i) can be 



from which exact values for g^"™(t) = 1 
determined. The resulting fitting constants Ocv, &cv and Cc 
are summarised in the top half of Table The method as- 
sumes the instantaneous removal of stars and ignores mass 
loss by stellar winds, but this assumption is legitimate since 
massive stars hardly contribute to overall cluster mass and 
low-mass stars only suffer significant mass loss during the 
last 10% of their lifetime. In order to determine /i'v^C^) and 
evaluate Eq. [5] also fJ-lv{t) is required, which will be dis- 
cussed separately in Sect. 12.2.21 

Since in the present study the stellar content of clusters 
is considered, we will be evolving the stellar mass func- 
tion of a cluster (and thereby indirectly the cluster mass) 
rather than evaluating cluster mass itself. Therefore, we 
will use mniax(i) from the isochrones instead of Eq. [5]to in- 
corporate stellar evolution. The time evolution of mmax(i) 
is shown in the left-hand panel of Fig. [T] for metallicities 
Z = {0.0004,0.004,0.02}. 



Mirit)+M!{it). 



(4) 2.2.2. Stellar remnants 



2.2.1. Luminous mass 

The first right-hand term of Eq. [3] can be expressed in 
terms of a mass fraction nl°^{t) = M'f*(t)/Mci,i that is 
still present in the cluster if stellar evolution were the only 
mass loss mechanisn:(f|: 



dMy_ 
dt 



c{t) 

'C7(0)' 



(5) 



where the fraction C(i)/C(0) is included to correct the mass 
loss for a possible changing mass function normalisation 
due to dissolution (which will be described in Sect. 12. 3p . 
Consequently, the mass fraction lost due to stellar evolution 
is defined by (7*°* (t) = 1 — {t) . The remaining total mass 
fraction iJ,l°*'{t) is the sum of the luminous mass fraction 
fi]^{t) and the stellar remnant (sr) mass fraction fJ.f^{t): 



(6) 

iLamers et al.l ()2005al ) have shown that for luminous cluster 
1 — fi]^"^{t) can be approximated as 



mass, g;';"'(t) 



logg, 



lum 
ev 



(logt 



(7) 



for log t > flcv , and where Oov , &ov and Ccv are constants de- 
termined by the IMF and metallicity of the cluster. We use 
the Padova 1999 isochrones (BerteUi et al. 1994, AGB treat- 
ment as in Girardi et al. 2000) to determine the maximum 
stellar mass TOmax(i) that is still present in the cluster at 
time t. This also provides the initial maximum stellar mass 
JTimaxi, which is taken to be the maximum mass at the 



^ As shown by 'Mie ske fc Baumgardtl l|2007f ) , the capture of 
stars by clusters is ineffective. 

^ If cluster mass also decreases due to dissolution, then the 
mass loss would be larger, see Sect. 12.31 



Stellar remnants can constitute a large fraction of the total 
cluster masf[E|. Especially during the later stages of cluster 
lifetime, the cluster content is likely to be dominated by 
remnants. As will be shown in Sect. [4] several observables 
can be significantly affected, making it essential to include 
remnants when studying cluster evolution. 

The distinction between total and luminous cluster mass 
can be made by writing M^'i°* = + M^^, where M^^ 

denotes the part of the total cluster mass constituted by 
stellar remnants. The evolution of dark cluster mass can be 
expressed as 



dM^ 
dt 



dM^ 

~dr 



dM^ 
dt 



(9) 



dis 



where at the right-hand side the first (positive) term de- 
scribes the production of remnants due to stellar evolution, 
and the second (negative) term represents remnant loss due 
to dissolution. Stellar evolution removes stars at the high- 
mass end of the mass function. Stellar remnant production 
can be included by assuming an initial-remnant mass rela- 
tion and leaving a remnant m ass upon such removal. 

From iKalirai et al.l (j2008l ). for white dwarfs (mg < 
8 Mq) this relation is given by 



= 0.109 + 0.394 Mr. 



(10) 



^ Values for the cas e of a Salpeter IMF are also pro- 
vided by ILamers et al.l l|2005al ) , but thes e are based on an 
older version of the GALEV mo dels (|Schulz et all l2002l : 
I Anders fc Fritze-v. AlvenslebenI [20031 ). 

^ If stellar mass-dependent cluster dissolution is omitted, after 
12 Gyr stellar remnants constitute about 30% of the total cluster 
mass. Within this fraction, the summed mass ratios of black 
holes, neutron stars and white dwarfs are 1:3.5:9. These ratios 
change when mass-dependent dissolution is included. 
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Luminous cluster 


mass (gir) 








Salpeter IMF 


Kroupa IMF 


z 




(2cv 


bcv 


Ccv 






Ccv 


0.0004 


68.5101 


6.83 


0.316 


-1.824 


6.84 


0.298 


-1.648 


0.0040 


69.8779 


6.80 


0.313 


-1.844 


6.80 


0.295 


-1.667 


0.0080 


71.6802 


6.76 


0.309 


-1.853 


6.77 


0.291 


-1.674 


0.0200 


68.1211 


6.70 


0.308 


-1.872 


6.71 


0.290 


-1.691 


0.0500 


49.4481 


6.63 


0.314 


-1.897 


6.65 


0.297 


-1.714 


Total cluster mass (qlf) 




Salpeter IMF 


Kroupa IMF 


Z 




dcv 


feov 


Ccv 


dcv 


bcv 


Ccv 


0.0004 


68.5101 


6.93 


0.271 


-1.855 


6.93 


0.255 


-1.682 


0.0040 


69.8779 


6.89 


0.271 


-1.872 


6.90 


0.256 


-1.696 


0.0080 


71.6802 


6.88 


0.265 


-1.877 


6.88 


0.250 


-1.701 


0.0200 


68.1211 


6.82 


0.263 


-1.893 


6.83 


0.248 


-1.716 


0.0500 


49.4481 


6.78 


0.263 


-1.908 


6.79 


0.249 


-1.731 



Table 1. Top: Fitting values for a, 

^lum 



bcv and Co 



q"^'^"^{t) using the Padova 1999 isochrones (see text) at five different metallicities for Salpeter (0.1 M© 



for the luminous fractional cluster mass decrease due to stellar evolution 

< TOs < mmax.i) 

and Kroupa (0.08 M© < mg < mmax.i) IMFs. The values of the initial maximum stellar mass rTimax.i correspond to 
the maximum masses at the youngest isochrones (logi — 6.6). The maximum difference between the fits and the exact 
9gv'"(0 curves is 3% after 19 Gyr. Bottom: Same as above, but for total fractional cluster mass decrease due to stellar 
evolution, accounting for stellar remnants. The maximum difference between the fits and the exact ql°^{t) curves is less 
than 10% after 19 Gyr. 
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Fig. 1. Left: Maximum stellar mass m,nax(^) as a function of age for Z = {0.0004,0.004,0.02} (dotted, dashed and 
solid lines, respectively). Right: Produced stellar remnant mass msr{t) as a function of age for Z = {0.0004,0.004,0.02} 
(dotted, dashed and solid lines, respectively). The peak corresponds to the transition from neutron star production to 
white dwarf production, i.e., the lifetime of an 8 Mq star. 



with rUsr the remnant mass and mg the stellar mass, where 
rris is equal to the initial stellar mass in our models since 
instantaneous death is assumed. For nis < 0.45 Mq, Eg.fTUl 
is not valid since the remnant mass then exceeds the pro- 
genitor mass. However, this does not cause any problems 
because it only occurs for ages much larger than the Hubble 
time. 

For neutron stars (8 M© < mg < 30 M©), the initia l- 
remnant mass relation is taken to be (|Nomoto et al.in"988f) : 



= 3.63636 x 10"^ (wg - 8 Mq) + 1.02 Mq. (11) 



Stars with masses iris > 30 Mq result in a black hole and 
since no definitive models for black hole formation exist, the 
corresponding remnant mass is assumed to be 

TOgr = 8 Mq. (12) 

This value is in agreement with dynamical masses obtained 
from X-ray binary observations (e.g., £asares 2006). When 
the above initial-remnant mass relation is linked to the 
maximum stellar mass »Timax(i), the remnant mass TOgr(i) 
that is produced at time t can be obtained. This varies 
with metallicity; for Z = {0.0004, 0.004, 0.02} examples are 
shown in the right-hand panel of Fig. [T] 

We can now compute the total cluster mass fraction in 
stellar remnants fJ-'^^(t) by integrating the remnant masses 
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rrisr over the mass function for the mass range of all stars 
that have ended their lives: 



(13) 



with Ns{0, nis) the initial form of Eq.[T] The total remaining 
mass fraction (t) is obtained by adding this to the mass 
fraction in luminous mass (see Eq.[6]). Similarly to (Zev™j 
provide fitting constants as in Eq. [7] for (7*°* = 1 — in 
the bottom half of Table [1] Again, we will not use these 
fits in the remainder of this study, since the stellar content 
of clusters is treated. Instead, the initial-remnant relations 
from Eqs.fTUl— [T^are used. Nonetheless, Table [1] is provided 
for convenience. 

Given the expression for /i^'^(t), the remnant production 
per unit time (first right-hand term of Eq. ^ can now be 
determined. By multiplying the time derivative of fil\,(t) 
with initial cluster mass and correcting the normalisation 
of the stellar mass function for a possible decrease due to 
dissolution, this fraction can be translated into the time 
derivative of non-luminous cluster mass 



dMg 



dt 



JWcl.i- 



'C(O)' 



(14) 



which is similar to the expression derived for the total clus- 
ter mass evolution due to stellar evolution (Eq. O. 

The above approach does not account for kick veloc- 
ities that are obtained by neutron stars and black holes 
upon their creation. The implications for integrated cluster 
properties are discussed in Sect. 16. 1[ point (3). 

2.3. Dissolution 

The effect of dissolution on the stellar mass function differs 
for clusters with or without the preferential loss of low-mass 
stars. In clusters that have reached energy equipartition, 
the low-mass objects tend to reside in the outer regions, 
while the most massive bodies are generally located near 
the cluster centre, and the cluster preferentially loses low- 
mass objects. In clusters without the preferential loss of 
low-mass stars, bodief0 of different masses are lost with 
almost equal probabilities. Both cases are treated in this 
section. 

The dynamical mass loss (second right-hand term of 
Eq. [3]) can be expressed as 



dM*°* 
dt 



dis 



Tdis 



1-7 



tnM, 



tot 7 



to 



(15) 



where Tdis is the dissolution timescale. In empirical stud- 
ies it was shown to be related to present clus t er ma ss as 
Tdis = to(Mci / M(7i)T b ylBoutloukos fc Lameri (|2003|) and 
iLamers et all (|2005a^ . iBoutloukos fc LamersI (l2003irfound 
7 = .62 ± 0.06, the same value as derived bv lGieles et al.l 
(I2OO4I) for the AT-body sim u lation s of clusters in tidal fields 
bv iBaumgardt fc Makinol (|2003D . Please note that A/*f' 
represents the total cluster mass and thus includes stars 
as well as their remnants. 

The rapidity of the exponential mass decrease in Eq. [15] 
is set by the dissolution timescale parameter to , which is es- 
sentially the dissolution timescale for a hypothetical 1 Mq 



cluster. Because the total disruption time iJi'?*''' of a cluster 
is determined by dynamical dissolution and stellar evolu- 
tion and thus depends on the adopted IMF and on metal- 
licity, we prefer the use of over ij^"*''' to characterise the 
strength of dissolution alone. Throughout this paper we will 
give the typical total disruption times associated with used 
values of the dissolution timescale to- 



2.3.1. Mass loss from clusters without the preferential loss of 
low-mass stars: the 'canonical mode' 

Cluster mass evolution can be computed numerically from 
its time derivatives due to stellar evolution and dissolution 
(Eqs. [5] and [T5)) . For a more detailed and more accurate 
description of cluster evolution we turn to its stellar con- 
tent. We distinguish between stellar mass-independent dis- 
solution ('canonical mode') and the preferential loss of low 
mass stars ('preferential mode'). When considering mass 
loss in the canonical mode, i.e., the evolution of the stel- 
lar mass function of a cluster without the preferential loss 
of low-mass stars, the influence of mass loss due to stel- 
lar evolution and dissolution is twofold. Stellar evolution 
decreases the maximum stellar mass r7imax(0' while disso- 
lution decreases the normalisation factor of the mass func- 
tion C{t) (see Eq. [T|) because stellar masses are randomly 
distributed throughout the cluster and thus all bodies have 
approximately similar probabilities of being ejected. 

The normalisation of the mass function C (t) is directly 
proportional to luminous cluster mass M^^"(i) because the 
latter is obtained by integrating the mass function (Eq. [1]) 
over stellar mass. We define a parameter /prof (i) to describe 
the fraction of mass loss occurring in preferential mode (see 
Sect. I2.3.2P and a fraction fsr{t) of the mass loss to occur 
in the form of stellar remnants. This implies that the factor 
1 — /pref (0 denotes the fraction of mass loss that occurs in 
the canonical mode, while the factor 1 — fsr{t) represents 
the fraction of mass loss that takes place in the form of 
luminous mass (i.e., not in stellar remnants). For mass loss 
in the canonical mode, i.e., without the preferential loss of 
low-mass stars, we have /pref (0 — 0. Then the time deriva- 
tive of luminous cluster mass from a cluster in the canonical 
mode, i.e., without the preferential loss of low-mass stars, 
is written as 



dt 



dis 



- [l-/pref(0][l-/sr(0] 

M^™dC 

~ c{t) dF' 



dAf*f 
dt 



dis 



(16) 



where C oc M'™ 



and thus dlnC/dt = dlnM^p/dt leads 



to the last equality, and the label 'can' indicates we are 
dealing with mass loss in the canonical mode. Substituting 
the total mass derivative from Eq. [15] yields after rearrang- 
ing 



dC 
"dt 



-T- = -[l-/prcf(i)][l-/sr(i)] 



A.f*f 



1-7 



toM, 



lum 



at). 



(17) 



'Bodies' denotes either stars or stellar remnants. 



In the canonical mode, the infiuence of mass loss on the 
stellar mass function is completely described by this ex- 
pression. If either the cluster only loses stars of the lowest 
masses (/pref (^) = 1) or all mass is lost in the form of stellar 
remnants {fsr{t) = 1), we have constant normalisation C{t) 
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and mass loss affects the cluster content in different ways 
(see below and Sect. 12.3."^ . 

To include the loss of stellar remnants, we specify the 
expression for the fraction of mass loss that occurs in the 
form of remnants /sr(t)- For mass loss in the canonical mode 
the spatial distribution of bodies is random. This implies 
that the fraction of the total mass lost due to dissolution 
that is lost in the form of remnants is equal to the ratio of 
remnant mass to total cluster mass 



(18) 



We can use this parameter to describe how the cluster mass 
in stellar remnants changes due to dissolution (the second 
right-hand term of Eq. [S]): 



~dr 



[l-/prcf(t)]/sr(t) 



dis 



dA-f^'°* 
dt 



dis 



M, 



-[l-/prcf(t)]/..(t)^ 

to 



tot 1-7 



(19) 



where we use the description of the total mass loss deriva- 
tive due to dissolution from Eq. [151 Again, the factor 
1 — /prcf(t) represents the fraction of mass loss that takes 
place in the canonical mode, while the factor fsr{t) denotes 
the fraction of mass loss that occurs in the form of stellar 
remnants. This expression completely describes the influ- 
ence of mass loss in the canonical mode on the total cluster 
mass in stellar remnants. 



2.3.2. Mass loss from clusters including the preferential loss 
of low-mass stars: the 'preferential mode' 

From A^-body sim u lation s, it is shown by 
iBaumgardt fc Makinol (|2003D how energy equiparti- 
tion affects the evolution of the stellar mass function. 
From their study it is evident that most clusters exhibit 
the preferential loss of low-mass stars, making it a very 
i mportant mechanism when considering cluster evolution. 
InjL amers et al.l (j2006() , analytical models for the evolution 
of the mass function are presented which are based on 
these simulations. In their models, dissolution no longer 
induces a uniform effect on the mass function by decreasing 
its normalisation. Instead, preferentially low-mass stars are 
removed. This can be approximated by a gradual increase 
in the lo wer mass limit of th e stars present in the cluster 
mniin(t) (|Lamers et al.l I2006D . By using this description, 
the slope of the mass function remains unchanged as the 
loss of low-mass stars is incorporated by increasing the 
lower mass limit. This is done in such a way that the 
mean stellar mass is always compar able to the mean stellar 
mass in the iV-body simulations bv IBaumgardt fc Makinol 
(|2003f ). in which the shape of the mass function does 
change due to the preferential loss of low-mass stars. 

To describe mass loss in the preferential mode, i.e., mass 
loss including the preferential loss of low-mass stars, we use 
the parameter /prcf (0 that represents the fraction of cluster 
mass loss that occurs in preferential mode at time t. The 
formulation allows for intermediate cases of energy equipar- 
tition, since /pret(t) can hav e any value between and 1 . 
The iV-body simulations by ' Baumgardt fc Makinol (j2003t ) 
suggest a rapid transition from a random ejection of bodies 
to the preferential loss of low-mass stars, resulting in almost 
a step function for /prof(t)- For a cluster that initially does 



not preferentially lose low-mass stars but reaches complete 
energy equipartition at t — tpiof we can then write 



/prcf(0 



for t < ipref 

1 for t > tpref- 



(20) 



However, because complete energy equipartition (/prcf (t) = 
1) is unlikely to occur (see Sect. S]), we will probably have 

fprci{t) < 1 for i > ipref. 

By explicitly integrating stellar mass over the mass 
function (Eq. [1]), the luminous cluster mass is obtained. 
Because in the preferential mode dissolution only affects the 
minimum stellar mass, we can write for the time derivative 
of luminous cluster mass: 



d7\fip 
dt 



prcf 



dis 



/prcf(t)[l-/r'W] 



dt 



dis 



C{t)r,{r 



2 - /?(TOmin) dt 



(21) 



for /3(mniin) 7^ 2 and 

d M^"™ \ "^'"^ / d M*°' 



dt 



dt 



-C(i)?7(m,„in)— In(rnmin), 



(22) 



if /3(minin) — 2, where P{mnun) is the slope of the stel- 
lar mass function (see Eq. [1]) and the label 'pref indi- 
cates we are dealing with mass loss in the preferential 
mode. The fraction of mass loss occurring in the preferential 
mode is represented by the factor fprcf{t), while the factor 
1 — /U.'''*(t) denotes the fraction of mass loss in the prefer- 
ential mode occurring in the form of luminous mass (i.e., 
not in the form of remnants, see below for more details). 
Again, we substitute the total mass derivative (Eq. [T5)) . 
which yields a differential equation describing the effect of 
the preferential loss of low-mass stars. After rearranging 
the terms, mass loss in the preferential mode thus results 
in an evolution of the lower mass limit m,„in(i): 



dTOn 



dt 



/p„,w[i-/.r'«)i-""'"^"'-'-">""""'"'" 



C{t)r]{m 



mill )io 



(23) 



for all values of /3(TOmin)- If a cluster is has not reached 
energy equipartition (/prof(t) = 0), mass loss is independent 
of stellar mass. In that case, the time derivative dm^in/dt — 
and we thus have constant minin(t)- 

We now describe stellar remnant loss from clusters that 
have reached energy equipartition. If a cluster is completely 
mass-segregated, remnants are produced in the cluster cen- 
tre and only become available for dissolution if remnants 
are the least massive bodies in the cluster. Therefore, the 
fraction of mass loss taking place in the form of stellar rem- 
nants that is used here f^^''°^{t) differs from the expression 
for canonical mass loss (Eq. [T5)l . 

From Fig. [1] we know that the produced remnant mass 
nearly always decreases with time, while m,„in(i) is a 
monotonously increasing function of t. This implies that 
there is a time tsr at which mmin(t) increases to a value 
larger than the smallest stellar remnant mass msr(ij3. For 



^ Incidentally, because the produced stellar remnant mass 
nearly always descreases with time, the smallest remnant mass 
at time t is the stellar remnant mass that is produced at that 
time. 
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all t > tsr mass can be lost in the form of remnantfl 
For these values of i, all remnant masses rnsr(0 — '^sr < 
'7imin(i) SkTB Considered to be immediately available for dis- 
solution (see Appendix [XI for a justification of this assump- 
tion). Consequently, they are immediately lost by dissolu- 
tion, because they are the least massive bodies in the clus- 
ter. 

The fraction of mass loss occurring in the form of stel- 
lar remnants is different for clusters with and without the 
preferential loss of low-mass stars. Therefore, we define a 
separate parameter /J^''''^(t) for the fraction of mass loss in 
the form of remnants if the cluster preferentially loses low- 
mass bodies. The time derivative of dark cluster mass (the 
second right-hand term in Eq. [5]) thus becomes: 

prcf 



dt 



dis 



dt 



dis 



= -/pref(i)/sr(i) 



prcf/.xf^ 



(24) 



which is nearly the same expression as Eq. [19] that is valid 
in the canonical mode. The fraction of mass loss occurring 
in the preferential mode is again represented by the factor 
/pref(i)- We have yet to specify the fraction of mass loss in 
the form of remnants for a the preferential loss of low-mass 
stars /P.'°^(t). If we consider a certain time interval dt, there 
are two possibilities: either the remnant mass available for 
dissolution (M^[''^'*', which is the total mass of all remnants 
with masses smaller than the lowest stellar mass present 
at time t) is so large that all mass loss during the interval 
dt can be accounted for by removing remnants only, or 
luminous mass has to be lost as well. In the former case, 
the minimum stellar mass f7imin(i) is not reached before 
the end of the time interval dt while losing the lowest mass 
bodies (i.e., remnants) and all mass loss takes place in the 
form of stellar remnants. For the latter case, TO„iin(0 
reached during the interval dt, and the fraction of mass 
loss to occur in the form of remnants is then described by 
the ratio of the available remnant mass and the total mass 
loss during the time interval dt. This implies 



1, 



/ dt 



(25) 



where the denominator of the second term between brackets 
can be obtained from the expression for cluster mass loss 
due to dissolution (Eq. llSp and the numerator is computed 
numerically from 



M, 



sr,dis 
cl 



(<) = / msr(ms)iVs(i, ms)dms 

"Imax(t) 



to 



dt' 



(26) 



Strictly spoken, this merely holds for (t) that intersect 
msi{t) only once. From Fig. [1] we see that at the transition from 
neutron star production to white dwarf production there is a 
possibility for curves of mmin(t) to increase to a value larger than 
msr(t), before briefly being overtaken again due to the change of 
produced remnant type. In that case, mmin(t) intersects msi{t) 
three times. However, the slope of mmin(t) is typically very steep 
compared to mar (4), and thus only a very small fraction of all 
cluster initial masses will pass through the described fluctuation 
over a negligible timespan. Therefore, we can indeed assume that 
for all t > tsT mass can be lost in the form of stellar remnants. 



The first right-hand term denotes the total produced rem- 
nant mass that is available for a given present mass func- 
tion Ns{t, nis) and increases with time. The upper integra- 
tion limit of the integral TOs™'"''''^(t) represents the initial 
stellar mass corresponding to remnants with mass equal to 
"^min (^ and the mass function is thus integrated for pro- 
duced remnant masses msr{t) < m^in(t). By introducing 
the second right-hand term in Eq. 1261 we subtract the part 
of the produced remnant mass that has already been lost by 
dissolution. This integral follows from the time derivative of 
non-luminous cluster mass (Eq. [M|) and describes all rem- 
nant mass that was lost in the preferential mode. Using the 
present mass function in the first right-hand term of Eq. [^Sj 
assumes that any change in the normalisation of the stellar 
mass function has a proportional effect on total remnant 
mass. The normalisation only changes for mass loss in the 
canonical mode (i.e., mass loss from clusters without the 
preferential loss of low-mass stars), for which dissolution 
is stellar mass-independent. Hence, if the normalisation of 
the stellar mass function were to change, the total mass in 
remnants would be affected accordingly, thereby justifying 
the above assumption. 

2.4. Total cluster evolution 

The description of cluster evolution from Sects. 12.11 to 12.31 
includes stellar evolution and four modes of mass loss: lu- 
minous mass loss and stellar remnant loss from clusters 
without the preferential loss of low-mass stars, and lumi- 
nous mass loss and stellar remnant loss from clusters that 
have reached energy equipartition and do preferentially lose 
low-mass stars. In this section the derived expressions are 
combined. 

The evolution of the stellar mass function is completely 
described by (1) the time derivative of its normalisation 
factor C(t) (Eq. [T7|) that describes the general decrease of 
the mass function, (2) the time derivative of the minimum 
stellar mass mnun{t) fEq. [23)) that describes the loss of low- 
mass stars, and (3) mynax{t) that describes the loss of high- 
mass stars due to stellar evolution. When combining these 
equations, the luminous cluster mass follows from 



it) 



At) 



ait) 



msNs{t, ms)dms 



(27) 



with Ns{t,ms) the stellar mass function from Eq. [T] 

Similarly, the expressions describing stellar remnant loss 
for clusters with and without the preferential loss of low- 
mass stars (Eqs. [T9l and [24]) can be combined with the ex- 
pression representing remnant production (Eq. [T4|) to write 



dMl _ dfj^ C{t)_ 
dt ~ dt ^'''C(0) 



(28) 



)i,j-totl-7 

- {[1 - /prcf(t)]/sr(i) + /pref(i)/r'W} , 

to 

with /^ov(^) from Eq. 1131 The first right-hand term de- 
notes the production of stellar remnants, while the sec- 
ond represents the loss of stellar remnants in the canon- 
ical mode (represented by the factor [1 — /prcf (0]/sr(i)) 
and in the preferential mode (represented by the factor 



^ These are the most massive remnants that are available for 
dissolution. 
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No preferential loss of low-mass stars With the preferential loss of low-mass stars 




7 8 9 10 11 7 8 9 10 11 

log Age (yr) log Age (yr) 

Fig. 2. Left Evolution of minimum and maximum masses of stars (solid) and stellar remnants (dotted) in a cluster 
with A/ci,i = 10^ M0, to = 3 Myr and without the preferential loss of low-mass stars. Right: Same as left, but including 
the preferential loss of low-mass stars. The tickmark at logt « 9.5 indicates the onset of the preferential mode (see 
Sect.SXT]). 



/prof(0/sr"'^(^))- This equation can be integrated for the 
total remnant mass in the cluster M^[(i). Finally, addition 
of the total remnant mass to the luminous cluster mass 
from Eq. [27] yields the total cluster mass as formulated in 
Eq.m 

The above set of equations can be solved numerically 
and represents a range of models that describe the complete 
evolution of cluster content for clusters with and without 
the preferential loss of low-mass stars. A simple recursive 
integration scheme is used. As a criterion for total clus- 
ter disruption, we use a lower luminous cluster mass limit 
M^™ = 100 Mq, though other values can be adopted if 
necessary. 

We now briefly illustrate the time evolution of the mass 
ranges of stars and remnants in a lO*' Mq cluster. These 
are shown in Fig. [2] for models with and without the pref- 
erential loss of low-mass stars. Solid lines describe the up- 
per and lower mass limits of the stars in the cluster, while 
the minimum and maximum masses of stellar remnants are 
represented by dotted curves. The left panel shows the evo- 
lution without the preferential loss of low-mass stars, while 
for the right panel it is included. 

Without the preferential loss of low-mass stars, the min- 
imum stellar mass is constant since mass loss occurs by 
removing stars of all masses and thus by decreasing the 
normalisation of the mass functions of stars and remnants. 
Stellar evolution causes the maximum stellar mass to de- 
crease. The maximum remnant mass remains constant at 
the maximum remnant mass that was produced in the clus- 
ter. On the other hand, the minimum remnant mass de- 
creases as remnants of lower masses are produced. Until 
total cluster disruption, bodies from a broad mass range 
can be retained. The cluster is completely disrupted when 
the normalisation constant C{t) approaches zero. 

When including the preferential loss of low-mass stars, 
the minimum stellar mass starts to increase as soon as mass 
loss occurs in the preferential mode. The maximum stellar 
and remnant masses exhibit the same behaviour as for clus- 



ters without the preferential loss of low-mass stars. Initially, 
the same holds for the minimum remnant mass. However, 
when the minimum remnant mass reaches the minimum 
stellar mass, the lowest mass bodies in the cluster are stars 
and remnants. This leads to a combined evolution of the 
lower mass limits of both. The cluster is completely dis- 
rupted when the stellar mass limits meet. 

3. Computation of photometric cluster evolution 

Cluster photometry is compute d from the P a dova 1999 
isochrones, that are described in iBertelli et al.l ([1994) and 
are b ased on spectral energy distributions from Kurucz 
(1992?) , but use a treatment of AGB stars as in lGirardi et al. 
(j2000|). The photometry computation is accomphshed by 
direct integration of luminosities over the stellar mass func- 
tion for a given age and initial cluster mass. This approach 
allows for greater flexibility when including the preferen- 
tial loss of low-mass stars, since the evolving mass func- 
tion is explicitly included in the computation. If existing 
SSP models had been adopted, this would not have been 
the case because such models only include fading by stellar 
evolution for a flxed mass function. 

For a cluster of arbitrary age t and initial mass M^^i, 
the stellar luminosities of the two isochrones at ages ti, 

i = {1,2} closest to t are integrated over the computed 
mass functions at these ages, with ti < t and t2 > t. This 
results in total cluster luminosities ici,A(ti, Afci,i) for pass- 
band A. These luminosities are then interpolated to obtain 
Lc\,\{t, Mci^i) and converted to absolute magnitudes. 

For existing SSP models, which only account for the 
effect of stellar evolution and therefore do not treat clus- 
ters near their total disruption, the above procedure suf- 
fices to determine photometric cluster evolution. However, 

ii ti < t < tdis < ^2 the cluster is disrupted before ^2 and 
there are no stars left at ^2- This leads to inadequate lumi- 
nosity computations if the above interpolation is used. In 
that case the stellar mass function of the cluster at age t 
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is adopted for both ages ti and t2 and the mass range is 
shifted to fit the appropriate maximum stellar masses at 
these ages. After integrating the resulting two mass func- 
tions, we obtain two luminosities Lci^\{ti, Md^i) for each 
passband. Interpolation then yields Lci^\{t, Mci,i). The de- 
scribed method does not accurately reproduce the lumi- 
nosity contribution of the lowest stellar masses because the 
mass function is shifted to fit to -Miiiax(ii, ^2)- However, 
this leaves the resulting magnitude almost completely un- 
affected since the high-mass end of the mass function com- 
pletely dominates cluster photometry. 

Contrary to some existing SSP models like GALEV, no 
gas emission is included in our photometric models. Line 
emission is only important for clusters that contain mas- 
sive stars that emit ionising photons. This implies that 
our photometry can be considered to be accurate for ages 
t k,8 Myr for solar metallicity and t > 20 Myr if Z = 0.0004 
(jAnders fc Fritze-v. Alvensleben|[20M) . 

4. Photometric properties of clusters 

In this section we apply the models described in Sects. [2]to[2] 
to investigate the effects of the preferential loss of low-mass 
stars, stellar remnants, IMF and metallicity on the mass, 
magnitude, colour and mass-to-light ratio evolution of clus- 
ters. We computed our models for cases with and without 
the preferential loss of low-mass stars, with and without 
stellar remnants, Salpeter and Kroupa IMFs, and metal- 
licities Z = {0.0004,0.004,0.008,0.02,0.05}. Moreover, we 
considered initial cluster masses Md^i in the range 10^ — 
10^ Mq and dissolution timescales to of 0.1—100 My£3- 
As a result, cluster evolution for total disruption times 
t^Xil'^^ > 10 Myr has been computed for cluster ages be- 
tween 10 Myr and 19 Gyr (the upper age limit of the stel- 
lar isochrones). Models for a broad range of parameters 
are publicly available in electronic form at the CDS and 
also at http://www.astro.uu.nl/~kruijs, while predic- 
tions for specific models can be made by the first author 
upon request. The most important results of our models 
are discussed in this section. 

4.1. The efFects of model components 

Accounting for the preferential loss of low-mass stars and 
including the mass of stellar remnants both have their ef- 
fects in the framework of our models. In this section these 
effects are considered. 

4.1.1. The effects of the preferential loss of low-mass stars 

The preferential loss of low-mass stars induced by en- 
ergy equipartition and possibly also mass segregation (sec 
Sect. [1]) can be expected to have a significant effect on the 

This is the typical dissolution timescale range (e.g., 
iLamers et aLll2005bD . corresponding to ti'if' ~ 10**— lO" yr for 
a 10* Mq cluster. It can be easily checked by considering that the 
total disruption time is of the order of ijgs''' ~ to(M'ci,i/M0)'^. 
Throughout this section to — 3 Myr will be used, which is the 
mean value of this range and is the typical timescale for clus- 
ters in the solar neighbourhood. This value is typical of tidally 
dissolving globular clusters on circular orbits at 3 kpc from the 
Galactic centre, or clusters on e ccentric (e = 0.7) orbits with an 
apogalactic distance of 8.5 kpc (jBaumgardt fc Makincill2003l ). 



magn itude and colour evolution of clusters (jLamers et al.l 
|2006!), but also on their mass and mass-to-light ratio. 
Complete energy equipartition (/prcf(^) = 1) implies that 
only bodies of the lowest masses are lost, an implication 
that does not seem likely to be in accordance with reality 
for two reasons. First of all, external perturbations are not 
strictly confined to the very outer radius of a cluster due 
to internal cluster dynamics. Therefore, dynamical mass 
loss is not confined to the very outer layer of the cluster, 
which results in the loss of stars or remnants with masses 
above mnun{t). Secondly, complet e energy equipartition 
may not be reached by a cluster ([Baumgardt &: Makinol 
2003), inducing only a partial preferential loss of low- 
mass stars. Therefore, we tuned the evolution of mniin(^) 
in our complete mod e ls to the A'^-body simulations by 
iBaumgardt &: Makinol (|2003D so that their mean stellar 
mass evolution is similar. This allows us to determine a 
step function form for /prof(i) that leads to a mean mass 
evolution that corresponds best to its counterpart in the 
A'^-body simulations. This analysis results in 

f Cf^ - / * < (r,Q-\ 

Jprcflij - I Q 4 fgj. ^ > ^^^^^^ l^yj 

with tprof = 0.2tj°g'''. A more gradual evolution would fol- 
low the simulations somewhat better , but a step function 
serves as a good approximation (see ILamers et al.l (|2006f) 
and Sect. jSj). 

The assumption of constant ^prcf /^dj's'*' = 0.2 serves as 
a typical timescale after which the preferential loss of low- 
mass stars can become important. It can be justified by 
considering the mass-dependences of tprof and t^^*,'^^. The 
total disruption time scales with cluster mass similarly to 
the dissolution timescale, i.e., t^^^l'^^ oc Af"-^^. On the other 
hand, tprcf can be expected to scale with the half-mass re- 
laxation time tih, i-e., iprcf oc M'^^^r^^^. If we adopt a mean 
mass-radius relation rh oc mO io±o o3 (|Larsenll2004bl ) . this 
leads to iprof oc A/"-^*^, which is comparable to the mass- 
dependence of the total disruption time. We also considered 
models with tpiof = ^rh, assuming the same mass-radius re- 
lation, which indeed yields values for iprof that are similar 
to 0.2t*°*''i. 

Furthermore, we do not consider primordial mass seg- 
regation (tprof = 0) for the model runs that are presented 
in this section. However, from other model runs where we 
did set ipiof = we know that the effects of the preferen- 
tial loss of low-mass stars on cluster observables are about 
10% stronger for primordial mass segregation than for the 
described results. Of course, tpicf can always be adapted to 
describe different forms of mass segregation. 

The formulation in Eq. implies that after tprof dy- 
namical mass loss simultaneously occurs in both modes. 
Hereafter, any reference to the 'preferential loss of low-mass 
stars' implies the use of this prescription for /pref in our 
models. The vali ie of .fmefit ^^pref) is different from the 
one presented bv ILamers et al.l (|2006 ^ because they did not 
compare its value to A^-body simulations: only tp^cf was 
treated in that study. 

In Fig. O the effect of the preferential loss of low-mass 
stars on the total cluster mass Af*[°*, F-band magnitude, 
V ~ I and mass-to-light ratio {M/Ly) evolution is shown 
for clusters with initial mass M^^i = 10^ M©, no inclu- 
sion of stellar remnants (i.e., all remnants are immediately 
removed), a dissolution timescale to = 3 Myr, metallicity 
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Age [Gyr] Age [Gyr] 

Fig. 3. Effect of the preferential loss of low-mass stars on (a) total cluster mass, (b) y-band magnitude, (c) V — I 
and (d) M/Ly evolution for clusters with initial mass Afd^i — 10^ M©, no stellar remnants, a dissolution timescale 
to = 3 Myr (t^°f ^ = 16—16.5 Gyr), metallicity Z = 0.02, and a Kroupa IMF. SoHd curves denote cluster evolution 
when the preferential loss of low-mass stars is included, while clusters without the preferential loss of low-mass stars are 
represented by dotted lines. The onset of the preferential mode tp^cf is denoted by a vertical line. 



Z = 0.02, and a Kroupa IMF. This leads to a total dis- 
ruption time of t^^^l^^ « 16 Gyr. For cluster evolution with- 
out the preferential loss of low-mass stars, some expected 
trends are immediately evident: mass and magnitude de- 
crease with time, while M/Ly increases and the cluster 
becomes redder. 

As can be observed when comparing the curves in Fig. [3l 
including the preferential loss of low-mass stars has several 
implications for the cluster evolution computed with our 
models. 

(1) The total disruption time of a cluster including the pref- 
erential loss of low-mass stars hardly changes but is 
slightly smaller than for a cluster for which it is omitted. 
This is the result of our model assumptions, as the pref- 
erential loss of low-mass bodies causes a larger number 
of massive stars to be retained in the cluster. The cor- 
responding shorter lifetimes induce an increase in the 
cluster mass loss by stellar evolution. The enhanced 
decrease of total cluster mass consequently leads to a 
smaller total disruption time. However, this effect does 
not include the possible decrease of the total disruption 
time due to quicker two-body relaxation in clusters with 
enhanced mean stellar masses. Note that the decrease 
of total disruption time is of order ~ 1%, which is best 
visible in the panel displaying colour evolution. 



(2) As soon as the preferential loss of low-mass stars starts, 
i.e., /prof(i) > 0, the cluster stays much brighter than 
for mass loss in the canonical mode. Because the lumi- 
nosity per unit mass is much higher for massive stars, 
a cluster that preferentially loses low-mass stars will be 
more luminous than a cluster of the same mass that does 
not. The change in the V^-band magnitude induced by 
the preferential loss of low-mass stars peaks at about 
0.9^^°*''' and is at most 1.5 mag. 

(3) The colour evolution of clusters is affected by the pref- 
erential loss of low-mass stars in two ways. After the 
onset of energy equipartition, Fig. \^c ) shows that these 
clusters are bluer than clusters without the preferential 
loss of low-mass stars, while just before total disruption 
reddening can be observed. At first, the bottom end of 
the main sequence, which is being lost due to the prefer- 
ential loss of low-mass stars, is redder than the average 
colour of the cluster. Due to the removal of its red con- 
stituents, such a cluster will appear bluer than a cluster 
without the preferential loss of low-mass stars. However, 
as stars are lost and m„iin(i) moves up the main se- 
quence the colours of the stars that are ejected become 
increasingly blue. Cluster colour is then dominated by 
red giants and the ejected stars are bluer than the aver- 
age cluster colour. Before total disruption, this induces 
a reddening of clusters exhibiting the preferential loss of 
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low-mass stars. The ages at which these changes occur 
depend on cluster tJi"*'*'. For model runs with different 
^dis'^' (either by varying initial cluster mass M^^i or dis- 
solution timescale Iq), we find that clusters including 
the preferential loss of low-mass stars always become 
slightly bluer from about OAt^^^l'^^ on. However, the red- 
dening is stronger and occurs at a smaller fraction of 
^dis'^' for clusters with smaller total disruption times. 
This tendency is caused by massive red giants, which 
are redder and more luminous than low-mass giants. 
For smaller total disruption times, the reddening causes 
clusters in the preferential mode to become much redder 
than clusters losing mass in the canonical mode during 
the last few percent of their lifetime. This effect is also 
present for the initial conditions used here if stellar rem- 
nants are included (see Fig. 2]). 
(4) As can be observed in Fig. [3Jd), the preferential loss of 
low-mass stars leads to a much smaller mass-to-light ra- 
tio than for clusters losing mass in the canonical mode. 
Considering points (1) and (2), this is not suprising. A 
higher luminosity and slightly smaller mass together im- 
ply a decrease in M/Ly. The magnitude of the induced 
decrease is comparable to the change in cluster F-band 
luminosity, which follows from the magnitude change to 
be about a factor six. 

In their study, iLamers et al.l (|2006[ ) obtained magnitude 
evolution curves that are much more weakly affected by the 
preferential loss of low-mass stars than is shown in Fig. [31 
In the present paper, cluster photometries are directly com- 
puted from the changing stellar mass func tion, which is 
a mor e direct method than the one used in iLamers et al.l 
(|2006l1^^l . Therefore, our computation of cluster evolution 
provides an update to their results. Furthermore, they 
found that clusters including the preferential loss of low- 
mass stars are bluer for 0.4 < < 0.8 and redder for 
t > 0.8tJ^°*^'. Our extended study of the parameter space 
shows that especially the value of 0.8 depends on the total 
disruption time of the cluster. 

The effect of the preferential loss of low-mass stars is 
clearly visible in magnitude, colour and M/Ly evolution 
of clusters. Especially cluster magnitude and M/Ly are 
significantly affected. This implies that when considering 
either observable, it is very important to detemine whether 
the cluster exhibits the preferential loss of low-mass stars or 
not, for example by checking whether it is mass-segregated. 
In order to obtain an appropriate interpretation of cluster 
colour, it has to be determined how close to total disruption 
the cluster is. 

4.1.2. The effects of stellar remnants 

The inclusion of stellar remnants follows the description 
presented in Sects. [^2.21 and [^31 In our models, it implies 
that part of the mass is retained upon the death of a star 
and clusters thus lose mass due to stellar evolution at a 
slower rate. It can be expected to affect cluster mass evo- 
lution, because a significant fraction of cluster mass can be 
constituted by remnants. Clearly, the mass-to-light ratio 
will then be altered as well. To assess the effect of stellar 
remnants on cluster evolution, our model results are shown 

ILamers et al.l (|2006l ) calculated magnitudes by using the 
photometry of stellar mass-truncated G^i£^F models. 



in Fig. [5] for clusters with initial mass Afci,i — 10^ M©, 
with and without the preferential loss of low-mass stars, a 
dissolution timescale = 3 Myr {t^^^'^^ = 16—18.5 Gyr), 
metalhcity Z = 0.02 and a Kroupa IMF. 

We observe the following changes in our models due to 
the inclusion of stellar remnants. 



(1) Because the net mass loss due to stellar evolution is 
smaller if stellar remnants are retained in the cluster 
after the death of their progenitors, the total cluster 
mass is higher than in the non-remnant case at any 
time. The immediate consequence is that the models 
predict a larger total disruption time. For clusters in- 
cluding the preferential mode, which more easily keep 
their remnants (see Sect. 12.3.21 and point (3)), this ef- 
fect is smaller than for clusters without the preferen- 
tial loss of low-mass stars, because the death criterion 
of clusters in our simulations (M'™ < 100 Mq) only 
depends on luminous cluster mass. Retaining remnants 
implies that mass loss due to dissolution more strongly 
affects luminous cluster mass, causing a cluster that ex- 
hibits the preferential loss of low-mass stars to satisfy 
the death criterion of our models earlier than expected 
for its total mass. However, though this weakens the 
lifetime-increasing effect of keeping stellar remnants, it 
never dominates. 

(2) For clusters without the preferential loss of low-mass 
stars, where the fraction of mass lost in the form of 
remnants is simply equal to the remnant mass fraction, 
the inclusion of remnants leads to a higher luminosity. 
Increasing luminosity by including remnants might be 
counter-intuitive. However, when including stellar rem- 
nants part of the cluster mass loss by dissolution is in the 
form of remnants instead of luminous stars. This implies 
that the average luminosity of bodies that are lost by 
dissolution is smaller than in the case without remnants, 
leading to a smaller luminosity decrease. Though less 
explicitly, the same effect is present in clusters includ- 
ing the preferential loss of low-mass stars since low- mass 
stars hardly contribute to the total cluster luminosity. 
Moreover, the lifetime-extending effect of remnants also 
implies that (luminous) mass is lost at a slower pace. 
We conclude that the cluster luminosity decrease due 
to dissolution becomes smaller if stellar remnants are 
included. 

(3) The mass-to-light ratio evolution shows a very clear ef- 
fect of stellar remnants as for both mass loss modes the 
M / Lv curves are much higher if remnants are included. 
This means that the relative increase of cluster mass due 
to remnants is larger than the corresponding relative in- 
crease of cluster magnitude that was discussed at point 
(2). The former is a direct consequence of adding rem- 
nants, while the latter is an induced effect: the average 
M/Lv of all bodies in a cluster is higher per definition 
if dark mass is added. Furthermore, for clusters with 
the preferential loss of low-mass stars that also includes 
stellar remnants, the mass-to-light ratio shows an in- 
crease during the final part of cluster evolution. This 
can be attributed to the preservation of remnants in 
clusters losing mass in the preferential mode, which is 
due to the fact that remnants can only be lost from 
these clusters if rus-c < TOmin(i)- 
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Fig. 4. Effect of stellar remnants on (a) total cluster mass, (b) F-band magnitude, (c) V ~ I and (d) M/Ly evolution 
for clusters with initial mass Md^i = 10^ Mq, a dissolution timescale to = 3 Myr {t^^^l^^ — 16 — 18.5 Gyr), metallicity 
Z = 0.02, and a Kroupa IMF. Solid curves denote the evolution for clusters including the preferential loss of low-mass 
stars with stellar remnants, long-dashed curves without stellar remnants. For clusters without the preferential loss of 
low-mass stars, short-dashed curves represent the case in which remnants are included and dotted curves represent the 
result without stellar remnants. The onset of the preferential mode tp^ci is marked by vertical lines in the linestyles of 
the corresponding model runs. 



Because remnant production or loss does not directly alter 
the colour composition of luminous cluster content, colour 
evolution is hardly affected by the inclusion of remnants. 

We find that the inclusion of stellar remnants strongly 
affects the mass, magnitude and M/Ly evolution of clus- 
ters. This is because part of the cluster mass loss occurs 
in the form of remnants rather than luminous stars. The 
extent of the differences (up to 30% at t ^ OM^^^l'^^ and 
increasing afterwards) suggests that a proper treatment of 
remnants should be included in any cluster model. 

4.2. The effects of the stellar IMF 

As discussed in Sect. 12.31 our models can be calculated for 
any multi-component power law IMF. Different IMFs are 
likely to exhibit a tendency to higher or lower stellar masses 
with respect to one another. For the Kroupa and Salpeter 
IMFs, the consequences of this effect are investigated here. 
As mentioned in Sect. [21 the lower mass limit is taken to be 
Wmin,i = 0.08 M0 for a Kroupa IMF and mmin.i = 0.1 Mq 
for a Salpeter IMF. Since a Salpeter IMF has a slightly 
steeper slope than a Kroupa IMF, and the latter features 
a bend at 0.5 Mq, the Salpeter IMF has a lower mean 
stellar mass. Figure [5] displays our results for both IMFs 
in the case of clusters with initial mass Afci.i — 10^ Mq, 



with and without the preferential loss of low-mass stars, 
with stellar remnants, a dissolution timescale to — 3 Myr 
{tf^^^ = 17.5—21 Gyr) and metaUicity Z = 0.02. 

The initial mass function affects the resulting cluster 
evolution in the following ways. 

(1) Total cluster mass stays higher for clusters with a 
Salpeter IMF. Consequently, their total disruption times 
are increased as well. Since it has a slightly steeper slope, 
the mass loss due to stellar evolution is smaller than for 
a Kroupa IMF. The resulting higher cluster mass leads 
to dissolution acting on a longer timescale, thereby also 
contributing to an extended lifetime of the cluster. In 
this case, the example model clusters with a Salpeter 
IMF survive beyond the maximum age spanned by the 
models. 

(2) Because it favours stars of higher masses, a Kroupa IMF 
leads to clusters that are initially slightly brighter than 
for a Salpeter IMF. However, the mass decrease due 
to stellar evolution is also stronger for a Kroupa IMF, 
causing its cluster mass to be lower than for a Salpeter 
IMF. As a result, the luminosities of clusters with a 
Salpeter IMF overtake those with a Kroupa IMF later 
on. Though still very weak, the transition is best visible 
when comparing the curves corresponding to clusters 
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Fig. 5. Effect of initial mass function on (a) total cluster mass, (b) y-band magnitude, (c) V — I and (d) M/Ly 
evolution for clusters with initial mass Mc\;i — 10^ M©, including stellar remnants, a dissolution timescale to = 3 Myr 
(^totai ^ i7_21 Gyr) and metallicity Z = 0.02. For a Kroupa IMF, solid curves denote the evolution for clusters including 
the preferential loss of low-mass stars, short-dashed lines for clusters with canonical mass loss. Clusters with a Salpeter 
IMF are represented by long-dashed lines when the preferential loss of low-mass stars is included and dotted lines describe 
the case where it is omitted. The onset of the preferential mode tprd is marked by vertical lines in the linestyles of the 
corresponding model runs. 



with mass loss in the preferential mode in Fig. [IJb). 
For these, the transition occurs at about O.Tt^'?*'*'. From 
model runs with other total disruption times, we observe 
that this fraction of tj°*'*' is not constant. It significantly 
mcreases for smaller t*^^l^\ with My of a Salpeter IMF 
always being fainter than that of a Kroupa IMF for 

^total < 5QQ ^yj.^ 

(3) Clusters with a Kroupa IMF are slightly bluer than those 
with a Salpeter IMF. However, this effect is too small to 
be observable in real clusters. It can be understood by 
considering the relatively larger contribution of massive 
(blue) stars in evolved clusters with a Kroupa IMF. Had 
the Salpeter examples been disrupted within the model 
age range, the characteristic reddening just before total 
disruption would have been visible for that IMF as well. 

(4) The higher masses of clusters with a Salpeter IMF and 
their generally lower luminosities that were described 
at point (2) together lead to higher mass-to-light ratios 
for clusters evolving from a Salpeter IMF. Again, for 
clusters including the preferential loss of low-mass stars, 
the final M/Ly increase upon total disruption would 
also be visible for the Salpeter IMF if it would have 
been completely disrupted within the model age range. 



Although specifically applied to the Kroupa and Salpeter 
IMFs, qualitatively the results of the above analysis hold 
for any two mass functions of which one has a different 
mean mass than the other. Quantitatively, there will still 
be variations depending on the specific IMF. 

Because the underlying stellar IMF in a cluster de- 
termines its future mass loss due to stellar evolution, it 
strongly affects the total disruption time and mass-to-light 
ratio. Any treatment of these two observables requires an 
accurate description of the IMF. 

4.3. The effects of metallicity 

To investigate the influence of metallicity on our results, we 
consider Z = {0.0004, 0.004, 0.02} in Fig. [H Model results 
are shown for clusters with initial mass Afci,i = 10^ Mq, 
with and without the preferential loss of low-mass stars, in- 
cluding stellar remnants, a dissolution timescale tg = 3 Myr 
(^totai ^ i7_i8.5 Gyr) and a Kroupa IMF. 

The effects of metallicity on the results are as follows. 

(1) We see that total cluster mass is hardly affected by 
metallicity at any time. It marginally increases with 
Z, which is caused by m ore rapid ste l lar evolutio n 
for low metallicities (e.g., iHurlev et all l2000l I2004D . 



14 



J. M. D. Kmijssen and H. J. G. L. M. Lamers: The photometric evolution of star clusters 




Age [Gyr] Age [Gyr] 

Fig. 6. Effect of metallicity on (a) total cluster mass, (b) F-band magnitude, (c) V — I and (d) M/Ly evolution for 
clusters with initial mass Mdj = 10^ Mq, including stellar remnants, a dissolution timescale to = 3 Myr {t^^l^^ — 17 — 
18.5 Gyr) and a Kroupa IMF. Solid curves denote the evolution for Z — 0.02, dashed curves for Z = 0.004, and dotted 
ones for Z = 0.0004. Results with and without the preferential loss of low-mass stars show effects as presented in Fig. |3] 
and described in Sect. I4.l7l1 and are therefore represented by the same linestyles. The onset of the preferential mode tprof 
is marked by vertical lines in the linestyles of the corresponding model runs. 



Consequently, the total disruption time also slightly in- 
creases with metallicity, which is best observed in the 
My and V — I panels. Nevertheless, the effect is small, 
less than 10% of the total lifeti me. This is in excell ent 
agreement with the results from lHurlev et al] ()2004[ ). 

(2) Low-Z clusters are brighter than high-Z clusters. This 
is due to a gen eral luminosity decrease of stars wit h 
metallicity (e.g. jGirardi et al.H2000HHurlev et al.l2004D . 
The difference is observed for clusters with and without 
the preferential loss of low-mass stars, and is typically 
more than one y-band magnitude. 

(3) The colour evolution shows a uniform trend with Z. 
Clusters with high metallicity are much redder than clus- 
ters with low metallicity. This is the res ult of stellar 
atmospheres and stellar evolution (e.g., iHurlev et al.l 
HOOO). For high ages, the V — I value is more or less 
constant for each metallicity, underlining its value as 
metallicity probe for glo bular clusters w hen considering 
broadband colours (e.g.. lMarastonl2005f ). The V—I shift 
between Z = 0.0004 and Z = 0.02 is about 0.5 mag, but 
varies for colours at other wavelengths. This is in accor- 
dance with the fact that clusters of different metallicities 
move on clearl y distinguishable pa ths in colour-colour 
diagrams (e.g.. lBruzual fc Charlotl f2003l. 

(4) The higher luminosity of low-metallicity stars and the 
consequently slightly enhanced mass loss by stellar evo- 



lution induce a common effect on cluster mass-to-light 
ratios. From Fig. EKd) we see that M/Ly strongly in- 
creases with metallicity. The effect is strong enough to 
move the M/Ly evolution of a Z = 0.0004 cluster with 
canonical mass loss through the M/Ly range of a clus- 
ter that does include the preferential loss of low-mass 
stars at Z = 0.02. This apparent degeneracy is lifted by 
taking cluster colours into account (see point (3)). 

Whenever cluster magnitude, colour and mass-to-light ratio 
evolution are considered, cluster metallicity plays an im- 
portant role. They are all strongly affected by the adopted 
metallicity. For globular cluster populations this implies 
that a range of colours and mass-to-light ratios can be cov- 
ered by a metallicity spread of the population or internal 
cluster processes like self-enrichment. 

In iLamers et al.l (|2005al) . an expression is provided for 
the total disruption time t^^^l'^^ as a function of ^o, Md^i and 
7. In this study, we find that depends on the inclusion 
of the preferential loss of low-mass stars, stellar remnants, 
IMF and metallicity. Therefore, values for tj"*''' are best 
obtained by integrating the models prese nted in t his pa- 
per. Regardless, the expression from Lame rs et al.l ()2005al ) 
can still be used to estimate i^"*''' with approximately 20% 
accuracy. 
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5. Application to globular clusters 

The results presented in Sect. 2] cover a wide range of 
masses, magnitudes, colours and mass-to-light ratios. As 
a first indication, it is relevant to check whether the prop- 
erties of Galactic globular clusters can be reproduced in 
our models. For this purpose, the results have t o be con- 
sidered at ages t w 12 Gyr. From iHarrisI (|1996D the My 
range is found to be My = -1.60 (Pal 1) to My = -10.29 
(wCen) . For the Solar neighbourhood value of the dissolu- 
tion timescale in the case of e = 0.7 orbits to — 3 MviF^. we 
find that this range can be covered at t = 12 Gyr for any 
metallicity Z < 0.02 if the maximum initial cluster mass 
equals M^'^^ = 10^ M©. For longer dissolution timescales 
(i.e., larger galactocentric radii) the observed magnitude 
range can be covered with even smaller maximum initial 
cluster masses. Please note that the dissolution timescale 
depends on the tidal field strength and that it therefore 
varies for different globular cluster orbits. This implies that 
it is not possible to impose limits on the properties of the 
complete globular cluster population (like their maximum 
initial masses) from an analysis in which this variation is 
not incorporated. 



5.1. The mass-to-light ratio 

Galactic globular cluster mass-to- light ratios are fou nd to 
be M/Lv = 1.45 ± 0.1 Mq Lq^ (iMcLauehlinl 120001) . SSP 
models, in which only stellar evolution is included and dy- 
namical effects are neglected, predict a minimum value of 
M / Lv ~ 2 Mq liQ^ at i = 12 Gyr, requiring the minimum 
metallicity of our models Z = 0.0004 (see Fig. This 
further complicates explaining the observed mass-to-light 
ratios, since gl obular cluster metallicities a re typically Z = 
0.0004—0.014 (jVanPalfsen fc Harrisl[200i) . Dynamical ef- 
fects are thus needed to explain the even smaller mass-to- 
light ratio of Galactic globular clusters. If the preferential 
loss of low-mass stars is included, cluster mass-to-light ra- 
tio curves do span the correct part of MjLy space, ranging 
down to MjLy < 1 Mq Lq^ (see also lKruiiss"enll2008D . 
From Fig. [5] we see that accounting for stellar remnants 
yields an increase up to 1 Mq Lq^ relative to model clus- 
ters without stellar remnants, obviously implying that they 
should be included for accurate interpretations of globular 
cluster observations. If globular clusters are populated us- 
ing a Salpeter IMF rather than a Kroupa IMF, this effect 
is nearly doubled (see Fig. [5]) . From Fig. [B]^d) we find that 
cluster models at high metallicities do not reproduce the 
low observed MjLy = 1.45 ± 0.1 Mq Lq\ Thus, metal- 
poor clusters including the preferential loss of low-mass 
stars with a Kroupa IMlQ are required to reach the low 
average mass-to-light ratio. If these conditions are met, the 
moderately flat MjLy evolution curves show that a more 
or less constant time-average is not surprising. 

Some globular clusters have mass-to-light ratios that 
are much higher than the mean value of MjLy ~ 1.45 ± 
0.1 Mq Lq\ An example i s wCen [MjLy = 2.5 Mq Lq^ 
see Ivan de Ven et al] |2006[ ) , which is also the most mas- 



sive Ga lactic globular clusteiH , with Af^^^:™ = 2.5 x 
10^ Mq ()van de Ven et al.1l200l . Since a high mass implies 
a large relaxation time, this could agrees with the view that 
wCen has not yet reached energy equipartition. 

Another important implication of our analysis of mass- 
to-light ratio evolution is that globular clusters of compa- 
rable ages can not be assumed to have constant MjLy 
for different cluster luminosities at fixed metallicity. For a 
given age, the mass-to-light ratio can strongly (~ 0.6 dex) 
depend on the dynamical state of the cluster (see Fig.[3I^d)), 
and thus on cluster mass and luminosity. The variation of 
MjLy with luminosity when including mass loss in the 
preferential mode is illustrated in Fig.[71 Models are shown 
of clusters with and without the preferential loss of low- 
mass stars, metalhcities Z = {0.0004, 0.004, 0.02} (dotted, 
dashed and solid curves, respectively), age t — 12 Gyr, a 
maximum cluster mass of 10 Mq, a Kroupa IMF, dissolu- 
tion timescale to = 3 Myr and including stellar remnants. 
The horizontal lines denote the constant mass-to-light ra- 
tios predicted if the preferential loss of low-mass stars is ig- 
nored, while the inclined curves show the relation between 
M/Ly and luminosity if it is accounted for. 

Since more luminous clusters are also more massive, the 
onset of mass loss in the preferential mode occurs later on 
for these clusters, implying that its effects are weaker for 
massive clusters at any age. Because the preferential loss 
of low-mass stars decreases the mass-to-light ratio, this de- 
crease is thus smaller for clusters of higher masses, leading 
to a mass-to- light ratio that increases with cluster mass and 
Ly as in Fig. [71 Observational evidence of this effect for the 
same quantitative range has been f ound for Galactic a,nd ex- 
tragalactic globular clusters (e.g., iMandushev et al.l 119911 : 



Rcikuba et aLll2007t lKruiissenll2008 '). If cluster masses are 
determined using a fixed M/Ly, thereby not accounting for 
the effects of the preferential loss of low-mass stars, these 
masses can be strongly overestimated by as much as 0.6 dex. 
Because the error is larger for lower masses, the slope of the 
inferred cluster mass function will be underestimated (i.e., 
a negative slope will be steeper) if the preferential loss of 
low-mass stars is ignored. 

It is straightforward to derive a quantitative estimate for 
the effect on the inferred (powerlaw) cluster mass function. 
Let us consider clusters with 'true' mass M that are exhibit- 
ing the preferential loss of low-mass stars and are thus in 
the regime that is inconsistent with canonical cluster mod- 
els. If we now use a powerlaw with index A to reasonably 
approximate the mass-to-light ratio increase with luminos- 
ity from Fig. [2l i.e., M/L cx L"^ for logL ^ 5.5 depending 
on metallicity, then the ratio of the photometrically inferred 
mass Mssp for which constant Mgsp/i is assumed to its 
true mass M scales as Mssp/Af cx L~^. This is equivalent 
to using a powerlaw with index B = A/(l -\- A) to approx- 
imate the mass-to-light ratio increase with true mass, i.e., 
M/L cx M^, leading to Mssp/M cx M"^. Then, for a 
'true' slope of the cluster mass function —a, its photomet- 
rically inferred slope — assp using constant mass-to-light 
ratios is given by 



assp 



-a-aA + A= {-a + B)/{1 - B). 



(30) 



12 See Sect. EXT] " in fact, wCen is not a normal globular cluster since there are 

1^ Or any other IMF that slightly favours massive stars with strong indications tha t it could be a stripped dwarf galaxy (e.g.. 



respect to a Salpeter IMF. Ildetafc Makinoll2004f ). 
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Fig. 7. Effect of the preferential loss of low-mass stars on the relation between mass-to-light ratio and luminosity at an 
age of t = 12 Gyr. Metallicities Z — {0.0004,0.004,0.02} are denoted by dotted, dashed and solid curves, respectively. 
Horizontal lines represent clusters without the preferential loss of low-mass stars, while inclined curves show the relation 
if it is included. The models are computed for initial masses between 10^ Mq and 10^ M©, a Kroupa IMF, dissolution 
timescale tg = 3 Myr and including stellar remnants. On each curve, clusters with (from left to right) initial masses 
logMci^i = {6, 6.5, 7} are marked with crosses (canonical mode) and dots (preferential mode). The present day mass can 
be derived from the luminosity and mass-to-light ratio. 



From Fig. [7] we find that typically A ~ 0.27 (and thus 
B ~ 0.21), implying that for a = 2 we find assp ^ 2.27. 
The deviation increases with the age of the cluster sample 
(here we used t = 12 Gyr, while for i = 3 Gyr assp ~ 2.23) 
and the above approach assumes a constant dissolution 
timescale for the cluster sample. For (more realistic) vary- 
ing dissolution timescales, the deviation will typically be 
between 50% and 100% of the the presented value. This 
is still significant and should thus be accounted for when 
studying the cluster mass function of old cluster samples. 

The variation of M/Ly with luminosity is crucial for 
cluster masses inferred from luminosities, which can be in- 
correctly determined by 0.6 dex. Because of increasing ev- 
idence that more globular clusters are preferentially losing 
low-mass stars and might have evol ved to a mass-segregate d 
state than previously thought (e.g. jDe Marchi et al.ll2007l ). 
this effect should always be considered when studying glob- 
ular clusters. 

5.2. Colours of globular dusters 

The typical colour range of Galactic glo bular clusters is 
0.75—1.1 in V"-/ (ISmith fc Stradeill2007D . From Fig.EIc) 
we can conclude that this range is covered by varying 
metallicity in our models. Obviously, IMF variations might 
slightly affect the colour range as well, with top-heavy 
IMF showing a tendency to bluer colours (see Fig. EI^c)). 



However, this effect is only significant for strongly deviat- 
ing IMFs. Still, from all figures in Sect. S] we still see that 
colours are affected by the preferential loss of low-mass stars 
and the colour range is increased. Figure [HIc) illustrates 
that if mass loss in the preferential mode would be ignored, 
uncertainties up to 0.2 mag (0.1 mag up and down) should 
be included for computed cluster colours at a any fixed age 
and metallicity. This happens t o be the same colour spr ead 
at fixed metallicity as found bv lSmith fc Straded (j2007t ) in 
their colour-metallicity relation. 

5.3. Ultra-compact dwarf galaxies 

We see that including the preferential loss of low-mass stars 
and stellar remnants allows for a more extensive analysis 
of globular clusters and affects their reproduced property 
ranges. It is evident that globular clusters or more mas- 
sive globular cluster-like objects with very high mass-to- 
light ratios (with regard to their metallicities) cannot have 
experienced significant mass loss in the preferential mode 
unless their MjLy is increased by a strongly differing IMF 
or by agents unaccounted for by our models, for instance 
by intermediate mass black holes, modified gravity or dark 
matter. Hence, if indications for the preferential loss of low- 
mass stars (such as mass segregation or a bottom-depleted 
mass function) are found for such objects, they points to 
these causes for the high observed M/ Ly- For objects with 
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masses M > 10^ M© this is unlikely to occur, since their 
relaxation times are of the order of a Hubble time or larger. 

If there is a present day globular cluster mass above 
which no clusters have yet reached energy equipartition, 
any objects above that mass will have high mass-to-light 
ratios with respect to clusters below that mass. This is 
found for UCDs, of which the M/Ly-range is typically 
2—10 Mq fe.g- lMieske fc Kroupal[2008l) . Pending the 
role of dark matter these galaxies could be regarded a nat- 
ural continuation of the mass spectrum beyond globular 
clusters. Indications for such a continuation are found by 
IWehner fc HarrisI (|2007t ) for UCD candidates in NGC 3311. 
As is clear from Fig. [Sj^d), our models without the prefer- 
ential loss of low-mass stars produce mass-to-light ratios 
similar to UCDs if a Salpetcr IMF is assumed (in agree- 
ment with Hilkcr ct al. 2007), while the canonical Kroupa 
IMF cannot reprod uce the range that is typical to UCDs 
(as is also found by iMieske et al.l [2008D . Nonetheless, any 
possible connection between globular clusters and UCDs 
would be smoothened by metallicity spreads and the pos- 
sible influence of dark matter, making it a challenge to be 
directly observed. 

Contrary to low-mass UCDs, more massive examples 
have not reached energy equipartition within a Hubble time 
due to their long relaxation times. However, this would 
lead to constant mas s-to-light ratios at high masses, while 
iReikuba et al.l (|2007t ) have shown that the trend of increas- 
ing M/Ly with mass is even stronger for UCDs than it is 
for globular clusters. 

6. Discussion 

We have described models to calculate the evolution of star 
clusters including the preferential loss of low-mass stars and 
stellar remnants. Furthermore, they have been used to in- 
vestigate the influence of these model components, as well 
as of the IMF and metallicity on cluster evolution. In this 
section we discuss the assumptions that were made, their 
influence on the results and the applicability of our models. 
We also indicate how the models can be improved. 

6.1. Influence of assumptions 

The models presented in this study are based on the fol- 
lowing assumptions. 

(1) We adopt the Padova 1999 stellar evolution models 
(Bertelli et al. 1994, AGB treatment as in Girardi et 
al. 2000) to compute photometry, determine stellar life- 
times and describe the consequent cluster mass loss due 
to stellar evolution. For the mass evolution of all stars 
we assume constant stellar masses until their instanta- 
neous deaths. Only very massive stars experience strong 
mass loss during a significant part of their lives, but 
these stars hardly contribute to the total cluster mass. 
Low-mass stars with mg ^ 8 M© only suffer significant 
mass loss during the last 10% of their lives. Therefore, 
instantaneous death is a reasonable assumption when 
calculating cluster mass. The photometric properties of 
stars are not affected by our assumption of instanta- 
neous death, since stellar photometry as described in the 
Padova models includes the photometric effects of stel- 
lar mass loss. The Padova isochrones include a descrip- 
tion of AGB evolution, which is very important since 



AGB stars dominate the photometric evolution of stars 
after about 100 Myr. 

(2) If stellar remnants are included, upon its death a star 
is replaced by a body with mass determined by the 
initial-rem nant mass rela t ion. F or white dwarfs, the re- 
lation from lKalirai et al] (|2008l ) is assumed, while neu- 
tro n star and bla c k hol e mas ses are b a sed o n studies 
by iNomoto et all (|1988D and ICasared (|2006[ ). respec- 
tively. Our resulting initial-remnant mass relation for 
the full stellar mass range is independent of metallicity. 
However, generally a me tallicity dependence is found for 
both white dwarf fe.g.. iKalirai et al ? 2005; Mcng ct ajj 
.2008. ) and neutron star masses (e.g.,lHurley et al, 200_0|). 
Remnants formed at high metallicities are generally less 
massive than those formed at low metallicities. The dif- 
ference between Z = 0.0004 and Z = 0.02 is typically 
^ 10%, which implies that the speed of the increase of 
non-luminous cluster mass M^f due to stellar remnant 
production would be enhanced by 10% at low metal- 
licities. For Fig. [HI which shows the effect of metallic- 
ity on our results, this would have some small conse- 
quences. Assuming a remnant to total cluster mass frac- 
tion M^[/M*°* « 0.5, the slope of total cluster mass evo- 
lution (Fig. El^a)) would be increased (i.e., become less 
steep) by a few percent at low metallicities, extending 
the total disruption time by a comparable but slightly 
lower percentage due to the M*°'''' mass dependence 
of the dissolution timescale. Moreover, curves describ- 
ing M/Lv cluster evolution (Fig. [Hd)) would exhibit 
metallicity effects that are smaller by a few percent. 
Overall, these corrections are not sufficiently large to 
have a significant effect on our results. 

(3) We assume all remnant bodies to be initially bound to 
the cluster. In reality, supernovae may induce kick ve- 
lociti es on their remnant black ho les and neutron stars 
(e.g.. IPortegies Zwart et anil997f ). A fraction of these 
remnants will have velocities larger than the escape 
velocity. Typical esca pe velocities of globular cluster s 
are ~ 30 km s^^ (M cLaughlin fc van der Marell [20051 ) . 
while in some cases kick velocities of several hun dreds 
km s~^ are observed (e.g. jLvne fc Lorimeilll99^ . This 
would imply that it would not be possible to retain 
all black holes and neutron stars in globular clusters, 
but nonetheless high numbers of neutron stars are ob - 
served in r eal globular clusters (jCamilo et al.l I2OOO0 . 
iPfahl et aD ([2b02il suggest that low kick velocities are 
obtained if neutron stars are formed in long-period and 
low-eccentricity high-mass X-ray binary (HMXB) sys- 
tems. In that case, the retained neutron star fraction 
would be four times higher than expected for commonly 
observed 'fast' neutron stars. Because of the lack of any 
definitive answer to the black hole and neutron star re- 
tention problem, and for the sake of model simplicity, 
we ignore kick velocities. In Fig. [8] we show the effect 
of including or excluding all black holes and neutron 
stars from our models. The panel is identical to that of 
Fig-Hl however curves representing model runs without 
stellar remnants are now replaced by ones which include 
white dwarfs only and for which all black holes and neu- 
tron stars have been removed. The effect is negligibly 
small (e.g., ^ 10% in all observables) because the total 
remnant mass is generally dominated by white dwarfs 
at ages ^ 100 Myr. In reality, the effect will be even 
smaller since a number of black holes and neutron stars 
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Fig. 8. Effect of neutron star ejection on (a) total cluster mass, (b) V^-band magnitude, (c) V — I and (d) M/Ly 
evolution for clusters with initial mass M^.i = 10^ M0, including all other stellar remnants, a dissolution timescale 
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represent the result where these remnants are removed upon their formation. The onset of the preferential mode t[pref is 
marked by vertical lines in the linestyles of the corresponding model runs. 



have velocities smaller than the escape velocity and are 
thus retained in the cluster. 

(4) Binaries are only partially incorporated in our mod- 
els. Dynamically, they are included since our mod- 
els are fitted to the collision al iV-body simulations by 
iBaumgardt fc Makinol (|2003 l). However, these do not in- 
clude a primordial binary population. Photometrically, 
binaries are not accounted for, because our model clus- 
ters are populated using single-star isochrones. As a re- 
sult, we have no mechanism in which white dwarfs can 
evolve towards neutron stars due to mass transfer. If 
we had, it would increase the total remnant mass by a 
very small amount proportional to the fraction of white 
dwarfs undergoing mass transfer in a binary system. 

(5) The preferential loss of low-mass stars is included by 
monotonously increasing the minimum stellar mass 
Wmin(i) of the bound stars in the cluster. As is dis- 
cussed in Sect. 14.1. li this is an approximation to the 
true evolution of the stellar mass function. T he iV-body 
simulations bv IBaumgardt fc Makinol ()2003D show that 
the slope of the IMF below a certain pivot-point mass 



mP'^ increasefl. thus exhibiting the preferential loss 
of low- mass stars. The maximum stellar mass ?nmax(^) 
is reduced by stellar evolution. Together, these effects 
narrow the mass function. In Fig. [51 our model mass 
and bolo metric mass-to-light ratio ev olution are com- 
pared to IBaumgardt fc Makinol (|2003D . The mass evo- 
lution shows good agreement, with a small offset at 
intermediate age that can be attributed to the dif- 
ferent stellar evolution prescriptions of both models 
(Pado va 1999 for the present paper a nd iHurlev et al.l 
(|2000l ) for lBaumgardt fc Makinol (|2003D ). Consequently, 
bo th models also differ in metall icity, as the value used 
by Baumgar dt fc Makinol (|2003l) is unavailable in the 
Padova isochrones. The difference in mass evolution is 
fully accounted for by these differences. 
The bolometric mass-to-light ratio evolution is very 
similar to our model including the preferential loss 
of low-mass stars for the largest part of cluster life- 
time. A strong difference between our approach and 
the iV-body simulations only arises in cluster photom- 
etry and remnant loss close to total cluster disruption. 



Because the slope is negative this means that it becomes less 
negative and eventually more and more positive. 
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Fig. 9. Comparison of the mass (left) and bolometric mass-to-light ratio evolution (right) from two of our models with 
and without the preferential loss of low-mass stars (Mci,i = 18407.6 Mq, with r emnants, tp = 22.5 Myr (t*,?^ ' ^ 11.5 Gyr), 
Z = 0.004 and a Kroupa IMF between 0.1 Moand 15 M©) to the results from lBaumgardt fc Makinol (|2003f ) (same initial 
mass, Wq = 5, Rgc = 8.5 kpc, circular orbit, Z = 0.001 and a Kroupa IMF between 0.1 MQand 15 Mq). 



when m^in{t) and m,max(i) nearly meet. In that case, 
model clusters only retain their giant branch, because 
increasing TOmin(i) first removes the main sequence be- 
fore giants are lost. This leads to luminosities that are 
overestimated close to total disruption. Moreover, this 
could cause real clusters to retain more remnants close 
to total disruption than is computed in our models. 
Consequently, the true mass-to-light ratio of clusters 
near total disruption would be larger than shown in 
this study. This is important in the last ~ 15% of total 
cluster lifetime (see Fig. [9]) , during which the mass-to- 
light ratio evolution can be expected to have a positive 
slope due to the retain of stellar remnants rather than 
the computed negative one. Therefore, the slight up- 
turn near total disruption induced by the inclusion of 
remnants that was shown in Fig. Sfd) can be expected 
to be much stronger in reality. Even a slight upturn 
does not appear in F i g. [H b ecause the IMF used by 
iBaumgar dt fc Makinol (|2003l ) is truncated at 15 Mq, 
implying that there is no production of black holes, 
which are needed for the upturn to occur in models 
where the preferential loss of low-mass bodies is rep- 
resented by an increasing lower mass limit. 
However, it is important to note that t — 0.85t|j°*'*' in 
Fig. [5] is also the point from where on the A^-body 
simulations are strongly affected by statistical noise. 
As the number of stars in the cluster decreases and it 
starts to be dominated by both high- and low-luminosity 
objects (giants and massive remnants) the imprint of 
statistics on especially the mass-to-light ratio is en- 
hanced. Consequently, it is very difficult to accurately 
describe the mass-to-light ratio evolution in the very last 
part of cluster lifetime, as such predictions are scatter- 
dominated and are thus not very likely to apply to any 
specific real cluster. 

The fractional range near total cluster disruption for 
which our model s do not follow the mass-to-light ratio 
rise from Baumgardt fc Makinol ([2003) increases with 
the total disruption time of a cluster, with 15% being 



the typical value for a t^^^^ of about a Hubble time. 
Shorter lived clusters are covered for by our models for 
a larger part of their lifetime, while for very long-lived 
clusters still at least 70% of their lifetime is covered. 
The vast majority of globular clusters is presently in 
the range where they can be treated with our models. 

(6) We adopt a step- function form for /pref(i) to describe 
the fraction of the cluster mass loss occurring in the 
form of low-mass stars if the cluster has (partially) 
reached energy equipartition. In the increasing mini- 
mum stellar mass approximation, a step-function re- 
quires the form of Eq. ^\ to most accurately repro- 
duce t he A^-body simulations by IBaumgardt fc Makinol 
(|2003f) . However, cluster magnitude evolution will be 
better reproduced if a smooth function of t is formu- 
lated. Nonetheless, this is not li kely to lead to an ex- 
act or better representation of IBaumgardt fc Makinol 
(j2003f ) due to the fundamentally different approach of 
including the preferential loss of low-mass stars. Instead, 
in a future study we will incorporate an improved de- 
scription of the changing mass function due to the loss 
of low-mass stars to account for the effect of the prefer- 
ential loss of low-mass stars (see Sect. 16. 2p . 

(7) Stellar remnants with masses smaller than the minimum 
stellar mass in the cluster, i.e., TOsi- < '7iiiii,i(t) are im- 
mediately available for dissolution. If the escape rate of 
bodies is not constant with cluster radius, this ignores 
the outward transport of remnants from their birth lo- 
cation (for a mass-segregated cluster this is considered 
to be in the cluster centre) to the cluster outskirts on 
the half-mass relaxation timescale irh- Because gener- 
ally < ^dis'^' — t except close to total cluster disrup- 
tion (see Appendix lA|l . this is a reasonable approach. 
In the exceptional case where ti± > t^^^l^^ — t, remnant 
loss is halted since the cluster is disrupted on a shorter 
timescale than trh- 

(8) The speed of cluster dissolution is assumed to be inde- 
pendent of cluster radius. GieLes_fc Baumgardt (200j) 
have shown that this is a reasonable assumption for 
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most tidally dissolving clusters, especially for large-A^ 
systems like globular clusters. 

6.2. Applicability and future studies 

Because this study is based on collisional A^-body simula- 
tions that confirm the existence of the preferential loss of 
low- mass stars and the retain of stellar remnants, the pre- 
dicted effects will be present in real clusters. It is shown 
that these phenomena, but also IMF and metallicity vari- 
ations can have unique effects on either the mass, magni- 
tude, colour and mass-to-light ratio evolution of clusters. 
Therefore, the effects of the preferential loss of low-mass 
stars, stellar remnants, IMF and metallicity can be ex- 
pected to be observable and interpretable. Clusters of all 
ages between 10 Myr and 19 Gyr can be treated with our 
models. Near total cluster disruption the results are affected 
by our formulation of the preferential loss of low-mass stars, 
in which red giants are the very last bodies to be lost near 
total cluster disruption. As a result, the cluster magnitude 
is overestimated and the mass-to-light ratio is underesti- 
mated during the last ~ 15% of cluster lifetime. Therefore, 
photometry-related observables have to be used with cau- 
tion near total cluster disruption (see Sect. 16. ip . 

The metallicity dependence of stellar remnant mass and 
an improved description of the preferential loss of low-mass 
stars are to be included in future studies. This will provide 
more accurate descriptions of non-luminous cluster mass, 
and cluster photometry and mass-to-light ratio close to to- 
tal disruption. A new set of evolutionary synthesis mod- 
els with a time-dependent stellar mass funct ion, based on 
A^-bo dy simulations will soon be available (|Anders et al.l 
|2008[) . Moreover, we have made a quantitative comparison 
of our models to a number of globular cluster systems, and 
have assessed u nexplained featu res of the mass-to-light ra- 
tio distribution (|Kruiissenll2008D . A paper in which the indi- 
vidual properties of Galactic globular clusters, like their or- 
bital parameters and metallicities, are used to explain their 
mass -to- light ratios is in preparation ((Kruiissen fc Mieska 
I2008D . 

7. Conclusions 

We have treated the influence of the preferential loss of 
low-mass stars, stellar IMF, metallicity and the inclusion 
of stellar remnants on cluster mass, magnitude, colour 
and mass-to-light ratio evolution. We presented analytical 
models that describe the evolution of cluster content and 
photometry, based on stellar evolution from the Padova 
1999 isochrones and on simp lified dynamical disso lution 
models as first presented in iLamers et al.l ()2005al ). The 
latter, in turn, is based on the A'^-body simulations by 
iBaumgardt fc Makinol (120031 ). 

The models represent the cluster evolution part of our 
new cluster population synthesis code SPACE. We con- 
sidered Kroupa and Salpeter IMFs and metallicities in 
the range Z = 0.0004 — 0.05. The obtained data are pub- 
licly available in electronic form at the CDS and also at 
http://www.astro.uu.nl/~kruijs. The results from our 
models are as follows. 

(1) The preferential loss of low-mass stars slightly decreases 
the total disruption time of a cluster by a few percent. 
However, the most significant changes are effected in 



cluster photometry. The effect of fading is decreased as 
clusters including mass loss in the preferential mode 
can stay more than 1.5 l^-band magnitudes brighter 
than clusters losing mass in the canonical mode, because 
most of the dynamical mass loss occurs in the form of 
low-mass stars that contribute little to cluster luminos- 
ity. Initially, clusters exhibiting the preferential loss of 
low-mass stars are bluer than standard ones, but they 
become redder during the last ~ 10% of cluster lifetime. 
The cluster mass-to-light ratio is severely decreased due 
to the preferential loss of low-mass stars. The decrease 
typically ranges from 2 — 4 M0 Lq^ (i.e., up to 0.6 dex) 
near total cluster disruption for total disruption times 
^totai ^ 12 Gyr. If the upturn o f the M/Ly evolution 
that is much more prominent in IBaumgardt fc Makinol 
(2003) than in our models (see Sect. 16.11 point (5)) is 
accounted for, this range of AI/Ly decrease is at most 
0.5 Mq Lq^ smaller. 

(2) Including the mass of .stellar remnants obviously yields 
an increase in the total cluster mass and consequently 
also in total disruption time with respect to cluster evo- 
lution without remnants. The extended lifespan also im- 
plies that cluster luminosity less rapidly decreases. The 
mass-to- light ratio is enhanced by almost 2 Mq Lq^ at 
its maximum, close to total disruption. 

(3) We compared the evolution of clusters with Salpeter 
and Kroupa IMFs, which can be considered to favour 
high stellar masses (Kroupa, or a 'top-oriented' IMF) 
or low masses (Salpeter, or a 'bottom-oriented' IMF) 
alternatives due to the bend in the Kroupa IMF and 
a slight slope difference. As can be expected, clusters 
with a bottom-oriented IMF retain more mass due to 
stellar evolution, which eventually causes these clusters 
to become brighter than clusters with a top-oriented 
IMF. However, they start out being slightly fainter 
since a top-oriented IMF favours massive stars and is 
thus brighter than a bottom-oriented one. Similarly, 
clusters with a top-oriented IMF are bluer and have 
smaller mass-to-light ratios than clusters with a bottom- 
oriented IMF. For the Kroupa and Salpeter IMFs, the 
latter change can amount up to several Mq Lq^. 

(4) Metallicity variations hardly influence the total mass 
evolution of clusters. In accordance with stellar stud- 
ies (Hur ley et al.l l2004f ) low-metallicity clusters are 
brighter and also much bluer than high-metallicity ones. 
Consequently, the mass-to-light ratio is a strongly in- 
creasing function of metallicity. 

(5) When applying our results to Galactic globular clusters, 
it is evident that the preferential loss of low-mass stars 
is required to explain their low observed mass-to-light 
ratios, especially if stellar remnants are accounted for. 
Low metallicity is insufficient to serve as an explanation. 
Another important implication of our study is that the 
mass-to-light ratio can not be assumed to be constant 
over varying luminosity, as it is strongly affected by the 
dynamical history of clusters. 

(6) The fact that clusters of high masses may not have 
reached energy equipartition yet suggests that the ef- 
fects of the preferential loss of low-mass stars disappear 
with increasing cluster mass. Because clusters exhibiting 
the preferential loss of low-mass stars have much lower 
mass-to-light ratios than clusters that lose their mass 
in the canonical mode, clusters with high masses would 
then have much higher mass-to-light ratios than ones 
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with lower ma s ses. T his effect may have been found by 
iReikuba et al.l (|2007f ). The above interpretation and its 
application to the observation s of.Reikub a et al. (20Q3) 



is treated more extensively in lKruijssenl ( 20081 ). 

(7) The typical colour range of globular clusters is covered 
by our models. When co nsidering the colo u r-met allicity 
relation as reported bv ISmith fc Stradeil (|2007| ). from 
an order-of-magnitude comparison we suggest that the 
observed colour scatter at fixed metallicity could be the 
effect of the preferential loss of low-mass stars. 

(8) Only when adopting a Salpeter IMF down to m^in,i — 
0.1 M0, the mass-to- light ratios of UCDs are reproduced 
by our models. While UCDs could represent a natural 
continuation of the trend of incre asing mass-to-light ra - 
tio with (globular) cluster mass (jReikuba et al.ll2007[) . 
this is not expected to be of a dynamical nature, since 
more massive UCDs are not expected to have reached 
energy equipartition within a Hubble time. 

The retain of remnants and the existence of the prefer- 
ential loss of low-mass stars are found in iV-body simula- 
tions of clusters and in observations, while metallicity and 
IMF variations are observed among real clusters. Therefore 
the effects described in this paper should be considered 
when observing clusters, and observed cluster properties 
have to be interpreted very carefullj{3- 
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Appendix A: Outward motion of stellar remnants 

In this study it is assumed that when stellar remnants are 
the lowest mass bodies in a cluster that preferentially loses 
low-mass stars they are immediately available for disso- 
lution. In mass-segregated clusters, remnants are created 
in the cluster centre, where the most massive stars reside. 
Because bodies are only lost from the cluster if they cross 
the tidal radius, this leads to a delay compared to remnants 
that would be produced at all radii unless the escape rate 
from the cl uster is independent of radius as proposed by 
iKind (|1966D . In this Appendix we show that the delay can 
be neglected even if the escape rate varies throughout the 
cluster. 

If the lowest mass bodies in the cluster are stellar 
remnants, these will move outwards on a half-mass re- 
laxat ion timescale Uh (e.g., ISpitzeH Il987t iHeggie fc HutI 
l2003l Ch. 14). Our description of remnant loss ignores any 
delay caused by the motion of remnants from the clus- 
ter centre to its outskirts. Therefore it implicitly assumes 
trh ^ ^dis'^' ~ where the latter term represents the re- 
maining lifetime of the cluster. To compare the two terms, 
we define % = (^1^°*'^' — If X > 1; remnants are able 

to reach the tidal radius before total cluster disruption; for 
X < 1 , the cluster is completely disrupted before such equi- 
llibrium can be reached. 

In Fig. lA.ll logx is shown for initial cluster masses 
between 10^ and 10'' Mq, with a dissolution timescale of 
to = 3 Myr, metallicity Z = 0.02, a Kroupa IMF and com- 
plete energy equipartition (/pref(^) = 1) after t = 0.2iJ^°*'^^ 
When a cluster reaches energy equipartition, there is a 
rapid drop in t^h because the mean mass in the cluster cen- 
tre increases. This is reflected in a sudden increase of logx 
with age for any specific initial mass. This can be observed 
in Fig. lA.ll at the dashed curve, which represents the onset 
of the preferential mode tprof (see Sect. 14.1. ip for the entire 
initial mass range. From Fig. lA.ll and model runs for other 
choices of dissolution timescale, metallicity and IMF, we 
can conclude that after iprof, X ^ 3 for all parameter sets 
relevant to globular clusters, and still x <^ 1-4 for clusters 
with initial masses Mdj < 10'^ Mq and extremely rapid 
cluster dissolution {to — 0.3 Myr). For massive clusters, 
we typically have x ~ 10 for t > ipref, implying that the 
immediate availability of remnants is a legitimate approxi- 
mation. 

Only during the very last few Myrs before total cluster 
disruption a cluster can have x < 1 because the numerator 
(remaining cluster lifetime) approaches zero more rapidly 
than the denominator (irh)- Remnants that are produced 
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Fig. A.l. Logarithm of the ratio of remaining cluster 
lifetime and the half-mass relaxation timescale, (t^"*'*' — 
t)/tj-ii = X, as a function of cluster age t and initial clus- 
ter mass Mci^i. The dashed line represents the onset of 
the preferential mode tp^ef for a cluster of corresponding 
y-axis initial mass, and the dotted line denotes the time 
from which on remnants are lost from such a cluster (tsr, 
see Sec. 12.3.21 for details). For the displayed model run we 
used a dissolution timescale of io = 3 Myr, metallicity 
Z = 0.02, a Kroupa IMF and complete energy equipar- 
tition for t > 0.2i^°f ^ 



during that brief phase cannot be lost in the preferential 
mode. If just before total disruption x < li the lower inte- 
gration limit minax(i) in Eq- [Ml that determines the total 
remnant mass available for dissolution, should be replaced 
by niinax(^x)' with t^ the time at which x decreases below 
unitjE!!- 

As we have shown, for other parts of cluster evolution 
the assumption of the immediate loss of remnants with 
masses msi(t) < mgr < m^init) is reasonable. Please note 
that for clusters including the preferential loss of low-mass 
stars with t > tprd, typically is not only very small 
compared to the remaining lifetime of the cluster, but also 
(^ 10%) compared to cluster age. 



To determine X; the total disruption time t^^l^^ has to be 
known before having integrated cluster properties over (and thus 
having obtained-) the entire cluster lifetime. This is alleviated 
by estimating t^^^^^ using straightforward integration of Eq. [3]at 
the very start of the computation. 



